The Riemann— Hilbert approach to strong asymptotics for 
orthogonal polynomials on [—1,1] 



A.B.J. Kuijlaars^ 

Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 200 B, 

3001 Leuven, Belgium 
arno@wis.kuleuven. ac.be 

K.T-R McLaughlin^ 

Department of Mathematics, University of North Carolina, Chapel Hill, NC 27599, 

U.S.A. 
and 

Department of Mathematics, University of Arizona, Tucson, AZ 85721, U.S.A. 

mcl@amath.unc.edu 

W. Van Assche^ 

Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 200 B, 

3001 Leuven, Belgium 
walter@wis.kuleuven.ac.be 

and 

M. Vanlessen^ 

Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 200 B, 

3001 Leuven, Belgium 
maarten.vanlessen@wis.kuleuven.ac.be 



Abstract 

We consider polynomials that are orthogonal on [— 1 , 1] with respect to a modified 
Jacobi weight (1 — + x)^h{x), with a,/3 > —1 and h real analytic and stricly 

positive on [—1,1]. We obtain full asymptotic expansions for the monic and or- 
thonormal polynomials outside the interval [—1,1], for the recurrence coefficients 
and for the leading coefficients of the orthonormal polynomials. We also deduce 
asymptotic behavior for the Hankel determinants and for the monic orthogonal 
polynomials on the interval [—1,1]. For the asymptotic analysis we use the steep- 
est descent technique for Riemann-Hilbert problems developed by Deift and Zhou, 
and applied to orthogonal polynomials on the real line by Deift, Kriecherbauer, 
McLaughlin, Venakides, and Zhou. In the steepest descent method we will use the 
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Szego function associated with the weight and for the local analysis around the 
endpoints ±1 we use Bessel functions of appropriate order, whereas Deift et al. use 
Airy functions. 



1 Introduction and statement of results 
1.1 Introduction 

Strong asymptotics for the classical orthogonal polynomials (Hermite, Laguerre and Ja- 
cobi polynomials) are known from the book of Szego [42, Section 8.21]. Szego himself 
developed the asymptotic theory of orthogonal polynomials on the unit circle, and on the 
unit interval for weights that satisfy the Szego condition, see [42] or [21]. During the last 
thirty years many asymptotic results were found for polynomials that are orthogonal with 
respect to various classes of weights, see the surveys [33, 34, 35, 39], and the monographs 
[31] and [44]. These developments were related to and partly motivated by questions 
from approximation theory, in particular Pade approximation and weighted polynomial 
approximation. 

Recent developments in the theory of random matrices [36] provided a renewed interest 
in orthogonal polynomials from a different perspective. For example, when the size of the 
matrices tends to infinity, local statistical properties of the eigenvalues are related to the 
strong asymptotics of the associated orthogonal polynomials. Important new techniques 
for obtaining asymptotics for orthogonal polynomials in all regions of the complex plane 
were subsequently developed by Bleher and Its [8] and by Deift et al. [10, 11, 12, 24]. The 
asymptotic analysis in the latter papers is based on the characterization of the orthogonal 
polynomials by means of a Riemann-Hilbert problem for 2x2 matrix valued functions 
due to Fokas, Its, and Kitaev [20], together with an application of the steepest descent 
method of Deift and Zhou, introduced in [15] and further developed in [4, 14, 16] and the 
papers cited before. See [10, 13, 25] for an introduction, and [5, 9, 17, 18, 19, 23, 26, 27, 
28, 29, 30, 37, 41, 45, 46] for the latest developments. A different method to obtain strong 
asymptotics for orthogonal polynomials was recently developed by Wong and co-authors, 
see e.g. [32, 40, 47]. 

The orthogonal polynomials considered in [11, 12] are orthogonal with respect to an 
exponential weight e^'^'^^^ on M or with respect to varying weights e~"^(^) on R. Inspired 
by these papers, we consider here polynomials that are orthogonal on a finite interval 
[— 1, 1] with respect to a modified Jacobi weight 



where a,f3 > —1 and h{x) is real analytic and strictly positive on [—1,1]. The main 
difference between the weights (1.1) and the weights e"'^^^^ and e~"^*^^^ where Q is a 
polynomial and V is real analytic on M, lies in the behavior near the fixed endpoints ±1. 
In the language of random matrix theory, this is the difference between the hard edge 
and the soft edge of the spectrum [38, 43]. While the analysis near the soft edge of the 
spectrum typically involves Airy functions, near the hard edge we expect Bessel functions. 

We use 7in{x) = 7r„(x; w) to denote the monic polynomial of degree n orthogonal with 
respect to the weight w on [— 1, 1]. Thus 



w{x) = {l-xy{l + xyh{x) 



X G [-1,1], 



(1.1) 




ior k — 0,1, . . . ,n — 1, 
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and Pn{x) = Pn{x;w) to denote the corresponding orthonormal polynomials. Thus 

Pn{x) = 7„7r„(a:;), 

where 7„ > is the leading coefficient of 

The weights (1.1) belong to the class of weights that satisfy the Szego condition 

^ \ogw(x) , ,^ ^, 

^ ^ ' dx> -oo. (1.2) 



-1 vr^ 



Szego 's theory, see [21, 42], provides the strong asymptotics of the polynomials Pn and 7r„ 
and the leading coefficients 7„ as n — > oo. To state the results, we need the function 

^^,) = , + ^,^-lYl\ ^eC\[-l,l], (1.3) 

which is the conformal map from C \ [—1, 1] onto the exterior of the unit circle. In (1.3) 
wc take that branch of {z^ — lY^"^ which is analytic in C \ [—1,1] and behaves like z as 
z ^ oo. Similar conventions will be used later on. We also need the so-called Szego 
function associated with a weight w satisfying (1.2), which for our purposes we define as 
the function D{z) — D{z; w) given by 

^/ X f(z^~iy^^ \0gw(x) dx \ , r 1 / .X 

i.W = exp(L_Ly_^_Sj^_j, f„r.eC\[-Ul. (1.4) 

The function D{z) is a non-zero analytic function on C \ [— 1, 1] such that 

D^{x)D_{x) = w{x), for a.e. x G (—1, 1), 

where D^{x) and D_{x) denote the limiting values of D{z) as z approaches x from above 
and below, respectively. The limit at infinity 

f 1 log'u;(a;) 
Doo = hm D{z) = exp I — / ax 



V27r y_i y/T^ 



X^ 



exists and is a positive real number. 

If the weight w on [—1,1] satisfies the Szego condition, then the following hold, see 
e.g. [42], [21]: 

li- S = ^ - 4 exp (-i /' . (1.5) 



n-oo 2" ^/7^D^ \ 271 J _^ ^/^ 



X 



n-.ooip{z)^ ^{z^-iy/W{z) ^ ' 



and 



hn, . ^^M!^ n 7) 

The limits (1.6) and (1.7) hold uniformly for z in compact subsets of C \ [— 1, 1]. 
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Remark 1.1 Since the Szego function of s/l — is equal to {z"^ — iy/^ip{z) we may 
also write (1.6) as 



lim 



1 



n^oo(^(^)" ^/2^D{z■,^/l-x^w{x)y 

which is the form usually found in the literature, see e.g. [44, Lemma 1.8.]. Because 
of (1.8), D{z] Vl — x^w{x)) is sometimes called the Szego function associated with w, 
instead of D{z;w). We found it more convenient to use (1.6). This form was also used, 
for example, in [6, 7]. 

Remcirk 1.2 There are also results for the orthogonal polynomials on the interval [—1, 1] 
under conditions that are somewhat stronger than the Szego condition, see e.g. [21, 31, 42]. 

In this paper, we give more precise asymptotic results than (1.5), (1.6), (1.7) for the 
special weights (1.1). We are able to obtain full asymptotic expansions for 7„, 7r„, and 
Pn, as well as for the coefficients a„ and 6„ in the three-term recurrence relation 



satisfied by the monic orthogonal polynomials. 

Prom our analysis we are also able to derive strong asymptotics for the orthogonal 
polynomials in the open interval (—1,1), as well as near the endpoints ±1. In a follow-up 
paper [29] we discuss its consequences for the level spacings of eigenvalues in the unitary 
random matrix ensembles associated with the weights (1.1). 

Besides the mapping function ip and the Szego function D, the statement of the results 
involves certain coefficients that are determined by the extra factor h in the weight. Since 
h is real analytic and strictly positive on [—1, 1], there is a neighborhood U of [—1, 1] such 
that h has an analytic extension to U (also denoted by h) with positive real part. Then 
log h is defined and analytic on U, where we take the branch of the logarithm which is 
real on [—1, 1]. For definiteness, we take U of the form 



for some < r < 1, where d{z, [—1, 1]) denotes the distance from z to [—1, 1]. 

Definition 1.3 Let 7 be a closed contour in U \ [—1,1], encircling the interval [—1,1] 
once in the positive direction. We define two sequences (c„) and {dn) of coefficients by 



TTn+i{z) = {z- bn)nn{z) - a^7r„_i(^). 



(1.9) 



U :^{zeC \ d{z, [-1,1]) < r} 



(1.10) 




(1.11) 



and 



]_ f log/i(C) 



dc 



(1.12) 



'n 



2m - 1) 



)V2 (C + 1) 



n+l ' 



forn = 0,1,2,.... 



Note that c„ and dn do not depend on the precise form of the contour 7. 
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1.2 Asymptotics for monic polynomials 



Our first result concerns the monic orthogonal polynomials 7r„. It is similar to Theorem 
8.21.9 of [42] where a full asymptotic expansion for the usual Jacobi polynomials Pi"'^\z) 
(that is, h = 1) was found. 

Theorem 1.4 For z E C \ [—1, 1], we have that 2'^nn{z)(f{z)~'^ has an asymptotic expan- 
sion in powers of 1/n of the form 



(p(zy ^ 



^{zy/^ 



D(z) y/2{z^ - 1)1/4 



k=l 



as n ^ oo. 



This expansion is valid uniformly on compact subsets of C \ [—1, 1]. The functions Ilk{z) 
are analytic on C \ [— 1, 1], and are explicitly computable. The first two are 



U,{z) 

Mz) 



4q;2 - 1 4/^2-1 
+ 



8((^(z) - 1) Si^piz) + 1) ' 
(4^2 _ 1) + /3 + Co) (4/32 -i)(a + /3 + do) (4a2 _ 1) (4/5^ - 1) 



1.13) 



+ 



16{^{z) - 1; 

2a'^ + 2(3^-5 
64 



16{^{z) + 1) 



128(^2 



4a' 



+ 



A/3' 



{^{z) - 1)2 {^{z) + 1) 



(1.14) 



Remark 1.5 The functions 11^(2;) are all rational function in ip{z), that vanish at 00. 
With increasing k, the functions Ilk{z) arc increasingly more complicated. 

If 4^2 — 1 = 4^2 _ _ then it turns out that all functions Ilk{z) are identically 
zero. In that case one can in fact prove that 



2"7r„(^) 



,1/2 



^{zy^ D{Z) ^2(^2 _ 1)1/4 

for some c > 0, cf. [25]. 



[1 + (e"'^")] , as n ^ 00 



1.3 Asymptotics for leading coefficients 

The leading coefficients 7^ of the orthonormal polynomials p„ also have an asymptotic 

expansion in powers of 1/n: 

Theorem 1.6 We have the following asymptotic expansion for 2~"'jn 



In 



1 



fe=i 



as n ^ 00, 



where the coefficients are explicitly computable. The first two are 
4q;2 - 1 4/32-1 



16 



16 



(a + /3 + Co) + — (a + /3 + do) 



32 



+ (2q;2 + 2/32 + 7) 



32 

4a2 - 1 4/32-1 
+ 



256 



256 



;i.i5) 



:i.i6) 



(1.17) 
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The coefficients are polynomial expressions in a, and the coefficients c„ and 
defined in (1.11) and (1.12). It turns out that and r2fc+i depend only on 

a, f3, Co, do, ... , Cfe_i, 

1.4 Asymptotics for orthonormal polynomials 

Since Pn — jnT^n, WB obtain from Theorem 1.6 and Theorem 1.4: 

Corollary 1.7 For z G C\[— 1, 1], we have thatpn{z)(p{z)~"' has an asymptotic expansion 
in powers of 1/n of the form 



Pn{z) 
(.p{zY 



1/2 



27r(^2 _ iy/^D{z) 



1 + 

fe=i 



as n ^ oo. 



The expansion is valid uniformly on compact subsets ofC\ [—1, 1]. The functions Pk{z) 
are analytic on C\ [—1,1], and are explicitly computable. The first two are 



Aa^ - 1 i^{z) + 1 _ 4/52 - 1 i^{z) - 1 
16 '^{z) - 1 16 (p{z) + 1' 



{4a^-l){a + (3 + Co)ip{z) + l , {Af]^ - l){a + P + do) ^{z) - 1 



+ 



32 

la^ - If 



+ 



ip{z) - 1 ^ 32 ip{z) + 1 

(4/^2-1)2 (4^2 -1) (4/52 -1) ip^{z) + l 



128{ip{z) - 1) 128((^(^) + 1) 

4n 2 - 1 



+ (2^2 + 2/3^ - 5) 
+ {2a^ + 2/3' + 7) 



64 

1 



6A{(f{z)-iy 6A{(f{z) + iy 



4^2 - 1 4/^2-1 
+ 



256 



256 



1.5 Asymptotics for Hankel determinants 

The Hankel determinant Dn associated with w is 

Dn ^det( f x^^^w{x)dx\ , 

\J-1 J j,k=0,...,n 

see [42, Chapter II]. Since 7„ = ■sj ^^'^ [42, Section 2.2], we have that 

n 



;i.i8) 



:i.i9) 



Prom Theorem 1.6 we then obtain the following asymptotic formula for Dn as n — > 00. 
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Corollciry 1.8 There is a constant C > such that 



as n ^ oo. 



Proof. Re- write (1.19) as 



-D„ _-i-r / 2^ \ 77/1,^11 TT I vTT^cx) 7j 



,rniv^o;:j ^^(l.^j ^Uy , 

The first two factors are expficit. The first one is equal to 



-2 



(1.21) 



)2 \ " 

oo 



2 J V 2 
and the second one is equal to 



(1.22) 



Cin-^ri (l + 0(i) ). (1.23) 



n 



The third factor is the partial product of a convergent infinite product, since by (1-15) 
we have 



2j 1 + ^i l + £i \f 



1 + 0^ as j ^ oo. 



If C2 is the value of the infinite product, then the third factor in (1.21) is 



Inserting (1.22), (1.23), and (1.24) into (1.21), and using (1.16), we obtain (1.20) with 
constant C = CiC2D{0). □ 

Remcirk 1.9 The 0{l/n) term in (1.23) and (1.24) can be expanded into a complete 
asymptotic expansion, with explicitly computable constants depending only on the coef- 
ficients Tj in the expansion of 7„. On the other hand, we cannot evaluate explicitly the 
constant C with our methods. For the determination of C in certain cases, see [22] and the 
references cited therein. This is known as Szego's strong limit for Hankel determinants. 

1.6 Asymptotics for recurrence coefficients 

The monic orthogonal polynomials satisfy a three term recurrence relation of the form 

TTn+liz) ^{Z- bn)TTn{z) - a^7r„_i(2;). 

In our next result we give an asymptotic expansion of these recurrence coefficients. 
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Theorem 1.10 We have the following asymptotic expansion for an and bn 
1 

k=2 



k=2 



as n ^ oo. 



The coefficients and are explicitly computable. The first few are 
4q;2 - 1 4/32 - 1 



32 



32 



(a + /3 + Co) + ^-^7^ — (a + /3 + do), 



32 

- (3q;^ + 3/3^ + 6aP + 1) 
4a2 _ 1 



32 

"4a2 - 1 4/32-1 
+ 



(1.25) 
(1.26) 

(1.27) 
(1.28) 



128 



128 



(4q;2_i)(4/32-1) 



256 



and 



Bo 



B. 



Ba = 



128 

13'' -a'' 
4 ' 

13'' -a' 
4 

Aa'-l 
64 



4/32 - 1 

3co(2a + 2/3 + Co) - 3do{2a + 2(3 + do), 



4q;2_1 4/32-1 
(1 + q; + /3) + Co — do 



16 



16 



(1.29) 

(1.30) 
(1.31) 



3co(2q; + 2/3 + 2 + Co) + 



4/32-1 
64 



3do(2Q; + 2/3 + 2 + 4) 



+ 



16 



(60; + 6/3 + 3q;2 + 3/32 + Gap + 4) . 



(1.32) 



Remark 1.11 Note that Ai = Bi = 0. The coefficients A2 and B2 only depend on a and 
/3. This means that up to order 1/^2 the recurrence coefficients agree with the recurrence 
coefficients of the pure Jacobi weight (1 — a;)"(l + x)^ asn— > 00. The effect of the extra 
factor h in w starts playing a role in the 1/n^-term by means of the coefficients Co and do- 
lt turns out that for the case 4^2 — l = 4/52 _ l = Q, all coefficients A^ and B^. vanish. 
In that case one can prove that a„ = | + 0(e~^") and 6„ = 0(e~'^") for some c > 0, cf. 
[25]. 

1.7 Asymptotics for monic polynomials on [—1,1] 

Our final results concern strong asymptotics for 7r„ on the interval of orthogonality [—1,1]. 
The asymptotic behavior of 7r„ in (—1, 1) is given by the following theorem. 
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Theorem 1.12 For x e (-1,1) 

TTnix) = 



2".0^(1 - x2)l/4 

^1 + 0(l/n)^ cos (^{n +1/2) arccos x + -0(2;) — 

+0(l/n)cos(^(n- l/2)arccosx + '0(a;) - ^) , (1.33) 

as n —> oo. The error terms hold uniformly for x in compact subsets of {—1,1) and have 
a full asymptotic expansion in terms ofl/n, which can he calculated explicitly. In (1.33) 
the phase function is given by 

,/ N 1/ / ^ n \ VI -x^ log/i(t) dt , 

-0(2;) = - Qi(arccosx — tt) + p arccosx H + , . (1.34) 

^ 2V ^ ' ^ ) 27r J_l^/^^^-x ^ ' 

The singular integral being understood in the sense of the principal value. 

The asymptotic behavior of TTn near the endpoint 1 is given by the following theorem. 
Theorem 1.13 There exists 6 > so that for x E (1 — 5, 1) 

^fixDao (narccosa;)^/^ 



7r„(x) 



2^*^^^ (l-x2)^/^ 

1 + 0(l/n)j ^cos^i(a;) Ja(narccosa;) + sin^i(a;) J^(narccosa;)^ 



X 



+ 0(l/n) (^cos C2(a;)<^a('T' arccosa;) + sinC2(a;) J^(narccosa;) j , (1.35) 

as n — > oo. The error terms hold uniformly for a; G (1 — 5, 1) and have a full asymptotic 
expansion in powers ofl/n, which can be calculated explicitly. In (1.35), is the usual 
Bessel function of order a, and 

1 CITT 

(1,2(2;) = ±- arccosx + ■0(a;) + — , (1.36) 

where the + holds for C,\, the — for (2, and where ip is given by (1.34). 

Near the endpoint —1, we can obtain similar asymptotic behavior in terms of Bessel 
functions of order /3. 

Remark 1.14 If a = ±1/2, it turns out, after an easy calculation using the facts 
Ji/2{z) = \ —sinz, J-i/2{z) = \ —cosz, 

' \l TTZ ' \l TTZ 

that the asymptotic formula (1.35) is the same as formula (1.33), so that (1.33) holds 
uniformly in [—1, 1] for these cases. This can be explained from the fact that for these 
special choices of a we do not need to construct a parametrix near the endpoint 1, see 
[25]. 
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Remark 1.15 The strong asymptotics (1.35) allow us to determine the asymptotics for 
the largest zeros of 7r„. Let 1 > Xn,i > Xn,2 > ■ ■ ■ denote the zeros of 7r„ numbered in 
decreasing order, and let < ja,i < ja,2 < • ■ denote the positive zeros of the Bessel 
function Jq, numbered in increasing order. Then, from (1.35) and Hurwitz' theorem we 
can deduce that, for every k — 1,2, . . . 



,■2 



= 1-^ + 0(1/71^), asn^oo. (1.37) 

This agrees with the well-known asymptotics for the largest zeros of Jacobi polynomials, 
see e.g. [2, Section 22.16] or [42, Section 8.1]. 

Remark 1.16 The strong asymptotic formula (1.35) is also related to the behavior of 
spacings of eigenvalues near 1 (the hard edge) of random matrices from the associated 
unitary ensemble, see [29] , where it was rigorously shown that the eigenvalue correlations 
can be expressed in terms of the Bessel kernel 

2{x - y) 



1.8 Overview of the rest of the paper 

As indicated in the beginning of the paper, our inspiration for this work comes from the 
papers [11, 12] by Deift, Kriecherbauer, McLaughlin, Venakides and Zhou. As in these 
papers, we use the Riemann-Hilbert approach to orthogonal polynomials, introduced by 
Fokas, Its and Kitaev [20] , together with the steepest descent method for Riemann-Hilbert 
problems of Deift and Zhou [15]. 

In this method a number of transformations Y ^ T ^ S ^ R are applied to the 
original Riemann-Hilbert problem in order to arrive at a Riemann-Hilbert problem for 
R, which is normalized at infinity, and whose jump matrices are uniformly close to the 
identity matrix. Then R is uniformly close to the identity matrix, which gives the leading 
terms in the asymptotic expansion for Y, and hence for the orthogonal polynomials. 
Using the fact that the jump matrices in the Riemann-Hilbert problem for R have a full 
asymptotic expansion in powers of 1/n, we arrive at a full asymptotic expansion for R. 
Tracing back the steps y h- > T i— > 5* h- > i?, we obtain the full asymptotic expansions for 
the orthogonal polynomials and related quantities given in Theorems 1.4, 1.6, 1.10, 1.12, 
amd 1.13 

All steps are analogous to the corresponding steps in [12]. The differences with [12] 
are connected with the endpoints ±1. The main differences are: 

• In the setup of the Riemann-Hilbert problem for Y we have to be careful with the 
endpoints ±1. IfQ;<0or/3<0, then the weight is unbounded at 1 or —1, and 
the Riemann-Hilbert problem cannot be stated in terms of continuous boundary 
values. We include a growth condition in the Riemann-Hilbert problem in order to 
control the possible growth of Y near the endpoints in these cases, see (2. 3) -(2. 4). 

• The first transformation Y t-^ T uses the conformal mapping cp from C \ [—1,1] 

onto the exterior of the unit circle. This step is simpler than the corresponding step 
in [12], since there the transformation is based on the construction of a (^-function, 
which depends on the weight and the degree n. There is also a preliminary step 
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Y \—>- U which scales the problem considered in [12], and which does not play a role 
here. 



• The second transformation T t-^ S involves a factorization of the jump matrix and 
a contour deformation. The factorization is somewhat different from the one in 
[12]. After contour deformation it leads to jump matrices for S with in some 
of the entries. These jump matrices therefore have a non-integrable singularity in 
case a > 1 OT /3 > 1. Consequently, the Riemann-Hilbert problem for S cannot 
be stated in L^-sense for some p > 1. Fortunately, the growth conditions at the 
endpoints still guarantee a unique solution, as we show. 

• Before we can do the third transformation S ^ Rwc have to construct a parametrix 
for S. The parametrix away from the endpoints depends in our case on the weight 
w. This is where the Szego function for w comes in. In [12] this particular step did 
not depend on w. 

• The construction of the parametrix near the endpoints ±1 involves (modified) Bessel 
functions (as expected), while in [12] it involves Airy functions. The construction 
is similar. One extra complication is that we have to take into account the growth 
condition at ±1. 

• Having the parametrix for S in all regions, we can do the third transformation 
S ^ R in the same way as in [12]. The asymptotic expansion of R follows from 
the known asymptotics of the modified Bessel functions at infinity, while of course 
it was based on the asymptotics of the Airy functions in [12]. 

Throughout the paper, w will be the weight (1.1) on [—1,1] with a,f3 > —1. The 
function h is positive on [—1,1], and analytic with positive real part in the domain 



with < r < 1. 

The function (f is the conformal map (1.3), D{z) is the Szego function (1.4) associated 
with w, and is the limit of D{z) for z ^ oo. Other notation will be introduced when 
needed. 

2 Riemann-Hilbert Problem 

The Riemann-Hilbert approach starts from a characterization of the monic orthogonal 
polynomial 7r„ with respect to a weight w on [—1,1] through the following Riemann- 
Hilbert problem (RHP) for a 2 x 2 matrix valued function Y{z). 

RHP for Y: 

(a) Y{z) is analytic for ;2 e C \ [— 1, 1]. 

(b) Y possesses continuous boundary values for x G (—1, 1) denoted by Y^{x) and 
Y„{x), where Y^{x) and Y_{x) denote the limiting values of Y[z') as z' approaches 
X from above and below, respectively, and 



U^{ze[-l,l]\d{z,[-l,l])<r} 




for X e (-1, 1). 



(2.1) 
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(c) Y{z) has the following asymptotic behavior at infinity: 



Y{z) 



I + 0{ - 

z 



z 





— n 



as z 



oo. 



(2.2) 



(d) Y{z) has the following behavior near z — 1: 



Y{z) 



O 
O 



1 Iz-l]"" 
1 Iz-V"' 

1 log 1-2—1 

1 log \ z — 1 



if a < 0, 
if q; = 0, 
if q; > 0, 



(2.3) 



^1 1 " 

as ^ ^ 1, ^ e C\ [-1,1]. 
(e) Y{z) has the following behavior near z — —1 



Y{z) 



^' 1 



o 



1 log 1 2; + 1 
1 log 1 2; + 1 

0,1 1- 



as z 



I 1," 
-1, ^eC\[-l,l]. 



if /3 < 0, 
if /3 = 0, 
if /? > 0, 



(2.4) 



Remcirk 2.1 The O-terms in (2.3) and (2.4) are to be taken entrywise. So for example 
Y{z) = C> (^J I Jjl) means that Yu{z) = 0(1), Fi2(;2) = 0{\z - 1|"), etc. 

Remark 2.2 The RHP for orthogonal polynomials is due to Fokas, Its, and Kitaev [20]. 
See Deift [10] for an excellent exposition. The above formulation of the RHP for weights 
on [—1, 1] differs from the formulation in [10] in two respects. 

• In [10] uniqueness of the solution of the RHP is formulated in the setting of 
boundary values. Here in (b) we assume continuous boundary values on the open 
interval (—1, 1), and we assume a certain growth condition at the endpoints ±1 in 
(d) and (e). If a > — | and > — |, then these growth conditions ensure that the 
RHP can also be taken in sense. This is no longer true ifQ;<— |or/3<— |. 
We will give an independent proof of the uniqueness of the solution of the RHP for 
Y in the first lemma. 

• The conditions (d) and (e) control the behavior at the endpoints ±1. There is no 
such condition in the RHP for orthogonal polynomials on the full real line, which 
are the object of research in [10, 11, 12]. We note that in the following we will 
encounter even more singular behavior at the endpoints. 



Lemma 2.3 If a solution of the RHP for Y exists then it is unique. 
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Proof. The proof is as in [10, p.44] except that we have to pay special attention to the 
behavior at the endpoints ±1. 

Let y be a solution. Then detY is analytic on C \ [—1, 1] and by the jump condition 
(2.1) we have for x e (—1, 1), 

{detY)+{x) = det (Y+ix)) = det (Y^x)) det (^^ ^^^^^ = (dety)_(x), 

so that det Y is analytic in C \ {—1, 1}. 
From (2.3) it follows that 

{0(|z-l|"), ifQ;<0, 
0(log|z-l|), ifa^O, 
0(1), ifQ;>0, 

as 2; — > 1. Since a > — 1, the singularity of detY at 1 is removable. Similarly it follows 
from (2.4) that the singularity at —1 is removable. Thus detY is an entire function. 

Prom (2.2) it follows that dety(2;) ^ 1, as 2; — > 00. Therefore, by Liouville's theorem, 
we have det Y{z) = 1 for every z. In particular, we see that Y[z) is invertible for every z 

and Y~^{z) is analytic for 2; e C \ [— 1, 1]. 

Now suppose that F is a second solution of the RHP for Y. Then H{z) := Y{z)Y^^{z) 
is defined and analytic for 2; G C \ [— 1, 1]. As in [10, p.44] we have that H is analytic 
across the interval (—1, 1). We have, since detY — 1, 

H(.) Y( AY-MA /^^ii(^) ^i2(^)\ ( Y22{z) -Y,2{z)\ 
H{z) = Y{z)Y (z) = ^^^^^^^ ^^^^^J ^^^^^^ j 

Yn{z)Y22{z) - Y,2{z)Y2i{z) Yi2{z)Yii{z) - Yn{z)Yi2{z) 
Y2i{z)Y22{z) - Y22{z)Y2i{z) Y22{z)Yn{z) - Y2i{z)Y,2{z) 

Since Y and Y both satisfy (2.3) it follows for each entry Hij of H that 

{0(1^-11°^), ifQ;<0, 
0(log|^-l|), ifa = 0, 
0(1), ifa>0, 

as z 1. Thus the point 1 is a removable singularity for each Hij. Similarly, the point 
— 1 is a removable singularity for each Hij. Therefore H is analytic in C Since H{z) — >■ /, 
as 2; 00, it follows again by Liouville's theorem that H = I and so Y = Y. 

This proves the uniqueness of the solution of the RHP for Y. □ 

Theorem 2.4 The matrix valued function Y{z) given by 
/ -n{z) ^J'^^-^dx \ 

V-27rz7Li^n-i(^) -7L1 ^-S^dxJ 

is the unique solution of the RHP for Y. 

Recall that 7r„ is the monic polynomial of degree n orthogonal with respect to the weight 
w and that 7„ is the leading coefficient of the corresponding orthonormal polynomial. 
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Figure 1: The contour 75 



Proof. It is obvious that Y{z) is analytic for 2; G C \ [—1, 1]. The proofs that Y satisfies 
(2.1) and (2.2) are as in [10]. We now concentrate on the proof of (2.3). It is clear that 
the (1, 1) and (2, 1) entries of Y are 0(1) as 2; ^ 1. For the (1, 2) entry, we have to prove 
that for z —>■ 1, 



0(|2-1|"), ifa<0, 
Yuiz) = { 0(log|z-l|), if « = 0, (2.6) 
0(1), ifa>0. 
This is obvious if a > 0, since then the limit 

\imYu{z) = — ^ ^ ' dx 2.7 
z^i 2m J_i x-l 

exists. Note that the integral in the right-hand side of (2.7) converges if a > 0. 

Consider the case a = 0. Then 

1 'Kn{x){l + xfh{x), , 1 /•i7r„(x)(l + a;)%(a:) , 

Yi2{z) = — / dx + — / dx. 

2m J_i X — z 2m Jq x — z 

The first term is analytic in a neighborhood of 1. For the second term, we note that 
7in{x){l + xYh{x) is analytic on [0,1], hence it satisfies the Holder condition on the 
interval [0, 1] so that by [1, p. 521] we have 

' ' ^^dx = c\og{z-l) + %{z), 



2m Jq X — z 

as 2; — i> 1, 2; G C \ [0,1], with some constant c and a bounded function $o(-2)- Thus 
Yi2{,z) = 0(log I2; — 1|) as 2 — 1, in case a = 0. 

Finally we consider the case a < 0. Recall that the function h is analytic in the open 
neighborhood U of [—1, 1], see (1.10) and Figure 1. Choose < p < r, so that the disk 
with radius p and center 1 is contained in U. Let \z — 1\ < p with z ^ [—1,1]. For 
6 > sufficiently small, z lies inside the closed contour 75 also shown in Figure 1. We 
define the function nn{C){C ~ 1)"(C + 1)^^(0 with a cut along {—00, 1], so it is analytic 
in U \ {—00,1]. Therefore, by Cauchy's theorem. 
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Letting 6 ^ 0, we find after an easy calculation that 

with Vt[z) analytic for \z — 1\ < p. Therefore 

sm ttq; sm ttq; 

Hence Yi2{z) — 0{\z — as 2; ^ 1, in case a < 0. 

So (2.6) holds in all cases. Similarly it follows that Y22{z) has the correct order as 
2; — > 1. Thus (2.3) holds. The proof that (2.4) holds follows in the same way. □ 

3 First transformation: Y ^ T 

In this section we formulate an equivalent RHP. This new RHP is normalized at infinity, 
and its jump matrix has rapidly oscillating entries. Let (T3 be the Pauli matrix 

A 
'"=^^(,0 -1 

and define 

T{z) = T"'^Y{z)ip{z)-'"'\ for ;2 e C \ [-1, 1], (3.1) 

with Y the unique solution of the RHP for given by (2.5). 

Recall that '^{z) = z + (z"^ — lY^"^ is the conformal map from C \ [—1, 1] onto the 
exterior of the unit disk. So it is clear from (3.1) that T is analytic in C \ [—1,1]. 

We have for x G (—1, 1), 

ip+{x)(p4x) ^ (x + Wl - a;2 j (^x - iVl - x^^ = 1, (3.2) 

and so, by (2.1) and (3.1), T satisfies the following jump relation on (—1, 1) 
T^x) = 2-'^^Y4x)^4x)--'^^ 

= 2-3r_(^)(^J '^f^)<^+(x)- 



-ncrs 



Using (2.2), (3.1), and the fact that <f{z) = 2z + 0(1/ z) as z ^ 00, we have T{z) = 
I + 0{l/z) as z —>■ 00. Finally, ip is bounded and bounded away from in a neighborhood 
of ±1. Thus, T has the same behavior as Y around the points ±1. So, we proved that T 
is the solution of the following RHP. 
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RHP for T: 

(a) T{z) is analytic for z e C \ [-1, 1]. 

(b) T{z) satisfies the following jump relation on 

w(x) 



T+{x) = r_(x) 







-2n 



-1,1): 

for X e (-1, 1). 



(c) T{z) has the following behavior at infinity: 



T{z) 



I + 0{- 

z 



as z 



oo. 



(3.3) 



(3.4) 



(d) T{z) has the following behavior as z 

O 



\z 
\z 



1|" 
11° 



o 



1 logl^; 
1 log|2; 

1 1 



O 



1 1 



^ 1: 

if q; < 0, 
if q; = 0, 
if a > 0. 

e) T{z) has the following behavior as 2; — > — 1: 

if /3 < 0, 

if /? = 0, 

if /3 > 0. 



(3.5) 



o 



1 + 

1 \z + l\^ 

1 log 1^; + 1 
1 log 1^ + 1 



(3.6) 



1 1 
1 1 



The RHPs for T and F are equivalent, since if T is any solution of the RHP for T 
then it similarly follows that 

y(;2) = 2-""8T(z)vp(^)"'"3^ (3.7) 

solves the RHP for Y. Thus, by Lemma 2.3, the RHP for T has a unique solution. 

Remcirk 3.1 Note that the transformation Y ^ T has the same form (3.1) for every 

weight w. What wc mean is that we use the same function ip for every w. This is in contrast 
to the case of varying weights on the real line [11, 12], where the first transformation is 
based on the construction of a so-called (^-function, which is specific for the weight in 
question (it may even be n-dependent). 
In our situation the g^-function is 

g{z)=\og^{z) = \og{z + {z' -l)"^), 

which is defined and analytic on C \ (—00, 1]. 

Remark 3.2 Since ^+{x) and ^-{x) have absolute value 1 for a; G (—1,1), we see in 
(3.3) that the diagonal entries in the jump matrix for T have absolute value 1, which are 
rapidly oscillating for large n. 
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Figure 2: The lens S 



4 Second transformation: T S 

The second transformation is based on a factorization of the jump matrix in (3.3) as a 
product of three matrices. A simple calculation, based on the fact that (p+{x)(p-{x) = 1 
for X G (—1, 1), shows that 

(p+{x)~'^"- w{x) 
(^_(a;)-2n 



1 0\ / w{x)\ f 1 

«;(x)-V-W"^" 1/ V-w(x)-i J lw;(x)-V+(a;)~^" 1 



(4.1) 



As in [10, 11, 12] we deform the RHP for T into a RHP with jumps on a lens shaped 
contour S = Si U S2 U S3 as shown in Figure 2. Here Si and S3 are the upper and lower 
lips of the lens, respectively, and S2 = [—1, 1]. The lens is contained in U. We also write 

S° = S,\{-1,1}, for J = 1,2, 3. 

The weight w has an analytic extension to U \ ((—00, —1] U [1, 00)), also denoted by 
w, given by 



w 



{z) = {i-zr{i+zfh{z). 



The extended weight is non-zero in f/ \ {{—00, —1] U [1, oo j), since Re/i > in f/. Let 
T be the unique solution of the RHP for T and define the matrix valued function S on 

C\S by 

T{z), for z outside the lens, 

S{z) = { T{z) (_^(^)-l^(^^)-2n, , for z in the upper part of the lens, (^4 2) 

T(z) ( / N 1^ / \ 9„ ? I , for z in the lower part of the lens. 

^ ^ \w{z) (p{z) ^" ly ' 

Straightforward calculations then show that S* is a solution of the following RHP. 
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RHP for 5: 

(a) S{z) is analytic for z e C \ E. 

(b) S{z) satisfies the following jump relations on S \ { — 1, 1}: 

S4z)^S4z)(,.,\,,_,„ °V forzeE?UE°, (4.3) 



w(;2)-V(^)"^" ly' 

S+ix)=S.ix){ ""[f^l, for a; e (-1,1). (4.4) 



-w{x)-^ 

(c) S{z) has the following behavior as z — > oo: 

S{z) ^ I + (^^^ , asz^oo. (4.5) 

(d) For a < 0, the matrix function S{z) has the following behavior as 2; — > 1: 

as ^ ^ 1, ^ e C \ E. (4.6) 





z — 


1 


(1 


z 


1 



For a — 0, S{z) has the following behavior as 2; — > 1: 

S{z) = oM'-]\ as.^l,.GC\E. (4.7) 

ylog|2; — 1| log|2; — l|y ^ ^ ^ 

For a > 0, S{z) has the following behavior as 2; — > 1: 

0||^ asz— >1 outside the lens, 
01 ^ , as2;— >1 inside the lens. 



(e) S{z) has the same behavior near —1 if we replace in (4.6), (4.7), (4.8) a by /3, \z — l\ 
by l-z + 1| and take the limit 2; — > — 1 instead of 2; — > 1. 

It is not completely obvious that the RHPs for T and S are equivalent. We have that 
the (unique) solution of the RHP for T gives rise to a solution of the RHP for S by means 
of the formulas (4.2). Conversely, we may start from a solution 5' of the RHP for 5" and 
define 

S{z), for 2; outside the lens, 

f(z) = { ^i^) (^w{z)-hp{z)-^'' 1) ' ^ "PP®^ P^^* (4-9) 

S{z) ^ ^i;(2;)~^(^(2;)~^" 1^ ' ^ lower part of the lens. 

We show that then T is a solution of the RHP for T, which in turn implies, since the 
solution of the RHP for T is unique, that the solution of the RHP for 5* is unique as well. 
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Lemma 4.1 Let S be a solution of the RHP for S and define T by the formulas (4.9). 
Then T is a solution of the RHP for T. 

Proof. It is straightforward to show that T satisfies parts (a), (b) and (c) of the RHP 
for T. We need to show that (d) and (e) also hold. Let T be the unique solution of the 
RHP for T. Then H{z) := T{z)T~^{z) is defined and analytic for 2; G C \ [— 1, 1]. Since 
T and T have the same jumps on (—1, 1), it follows that H is analytic across (—1, 1). We 
have, since detT = 1, 

ffu{z)T22{z)-f,2{z)T2,{z) fu{z)n,{z) -fu{z)T,2{z)\ 

H{z) = ( ^ ^ ^ . (4.10) 

\T2i{z)T22{z) - T22{z)T2i{z) T22{z)Tn{z) - T2i{z)T,2{z) J 

From the definition (4.9) of T{z) and (4.6)-(4.8) it follows that 

as z —>■ 1, 

\z — ■ 
in case a < 0, that 

^' Vlogl^; — 1| logl^; — l|y 





\z — 


ll°\ 


G 


\z — 





in case a = 0, and that 

f{z) = 










G 







- 1| 


—a 


(1; 


- 1 


—a 



as z ^ 1 outside the lens, 



O ( r "!_^ I 1 , as 2; ^ 1 inside the lens. 



in case a > 0. 

Combining this with (3.5) and (4.10), we obtain that 





\z — 


1 


a 1 


\z 


- 1! 




(1 


\z 


1 


a 1 


\z 


- 1 


F) 




as ^ ^ 1, (4.11) 



log I2; — 1| (log \z — 1\Y 
log l-s — 1| (log \z — l\Y 



1 1 
1 1 



asz^l, (4.12) 





\z — 




—a 


\z — 






(1 


\z 


1 


—a 


\z 


1 


1::) 



as 2; — > 1 outside the lens, 
as 2; ^ 1 inside the lens. 



(4.13) 



in case a > 0. 

In all cases it now follows that H has a removable singularity at 1. This is clear 
if a < from (4.11) and (4.12), since a > —1. If a > 0, then it follows from (4.13) 
that H remains bounded if we approach 1 from outside the lens. Then the entries of H 
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cannot have a pole at 1, since in the case of a pole, the function becomes unbounded if 
we approach 1 from all directions. Moreover, we also get from (4.13) that {z — l)"^H{z) 
is bounded near 1, for any integer m > a. Then the entries of H cannot have an essential 
singularity either, so that 1 is a removable singularity of H. 

Similarly, it follows that in all cases, H has a removable singularity at —1. Therefore, 
H is entire. It follows from (3.4), (4.5), and (4.9) that H{z) I as z ^ oo. So, by 
Liouville's theorem, H{z) = I. This gives T = T, so that indeed T is the (unique) solution 
of the RHP for T. □ 

Remark 4.2 For a > 1, the jump matrices (4.3) and (4.4) for S, as well as S itself, have 
a non-integrable singularity near z = 1, since the factor w~^{z) behaves like (1 — z)~°' as 
^ — > 1. It then follows that the RHP for S does not have a solution in any sense with 
p > 1, if a > 1. 

Remark 4.3 Since \(p{z)\ > 1 for every 2; e C \ [—1, 1], we see from (4.3) that the jump 
matrix for 5* on U S3, 

1 

w{z)-^cp{z)-^'' 1 

tends to the identity matrix as n ^ 00. So the effect of the transformation T ^ S has 
been to transform the oscillating entries on the diagonal in the jump matrix for T into 
exponentially decaying off-diagonal entries in the jump matrix for S. This lies at the 
heart of the Deift-Zhou steepest descent method for Riemann-Hilbert problems. 

5 The parametrix for the outside region 

Since the jump matrix in (4.3) tends to the identity matrix as n — > 00, we expect that 
the leading order asymptotics are determined by a solution N of the following RHP. 

RHP for N: 

(a) N{z) is analytic for 2; e C \ [-1, 1]. 

(b) N{z) satisfies the following jump relation on the interval (—1, 1): 

N^{x)=N_{x)( ° ''fA, for (-1,1). (5.1) 



-w{x) ^ 
(c) N{z) has the following behavior as 2; — >• 00: 

N{z) ^ I + (^^^ , asz^oo. (5.2) 

We present a solution of the RHP for N in terms of the Szego function D{z) = D{z] w) 
associated with w, see (1.4). The Szego function is the scalar function : C \ [— 1, 1] — > C 
that satisfies the following scalar RHP. 
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RHP for D: 

(a) D{z) has no zeros and is analytic for 2; e C \ [—1, 1]. 

(b) D+{x)D4x) = iLix) for x e (-1, 1). 

(c) lim D{z) = D^e (0,00). 

z—*oo 

It is well-known and easy to check that D as given by (1.4) is indeed a solution of this 
RHP. 

Remark 5.1 li Di and D2 are the Szego functions for the weights Wi and u»2, respectively, 

then the product -D1-D2 is the Szego function for W1W2. Since the Szego function for the 
pure Jacobi weight Wa,/3{x) = (1 — a;)"(l + x)'^ is equal to 

= for.GC\[-l,l], (5.3) 

it follows that for z e C \ [-1, 1], 

D{z]w) = D{z]Wa,i3)Diz;h) 

(^- l)"/2(2+ l)/3/2 /(^2_i)l/2 rl logh(x) dx \ 



^(^zY^+m 27r J.iVT^z-x 

Proposition 5.2 A solution of the RHP for N is given by 

/ a{z)+a{z)-'^ a{z)-a{z)-'^ \ 
V -2i 2 / 

where D{z) — D{z;w) is the Szego function associated with w, and 

= l^-p^. (5.6) 

Proof. If we seek a solution N of the RHP for N of the form 

N{z) = D2N^'^\z)D{z)-''\ (5.7) 
then Ar(i)(^) is analytic in C \ [-1, 1], and 

Ar«(^) = / + , asz^oo. 
Furthermore, using (5.1) and D^{x)D_{x) — w{x) for x e (—1, 1), we have 

N^+\x) = N^}\x) {^^^ , for X e (-1, 1). 
It is well-known (see [10, 11, 12]) that the RHP for iV(^) is solved by 



-1 / a{z)+a{z) ^ a{z)—a{z) 



-1 ■ 



2 2i 
a{z)—a{z)~^ a{z)+a{z)~^ 
2 

where a{z) is given by (5.6). This proves (5.5). □ 
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6 Parametrix near the endpoints 



The jump matrices for S and N are not uniformly close to each other near the endpoints 
±1. This is reflected in the fact that SN~^ is not bounded near the endpoints. We need a 
local analysis around these points. We consider the endpoint z = 1 in detail, the analysis 
for z = —1 being similar. In a small but fixed neighborhood Us = {z \ \z — 1\ < 6} of 
z = 1, with < 5 < r (so that the closure of Ug lies in t/) we construct a parametrix P, 
that satisfies exactly the same jump relations as S, that matches with N on the boundary 
dUs of Us, and that has the same behavior as S{z) near z — 1. We thus seek a matrix 
valued function P that satisfies the following RHP. 



RHP for P: 

(a) P(^) is defined and analytic for 2 e C/^g \ E for some So > S. 

(b) P{z) satisfies the following jump relations on Us n E°: 



P+(z) = P_(z) (^^^^^7^^^^-2n °), forzeC/,n(E°UE^), (6.1) 

P+(x) = P_(a;)(^_^°^^_, "^^^^y iorxeUsHEo. (6.2) 

(c) On dUs we have, as n — > 00 

P{z)N-\z) ^I + O , uniformly for zedUs\ E. (6.3) 

(d) For a < 0, the matrix function P{z) has the following behavior as 2; — > 1: 



a 



P{z)^0{] \l Jl ), asz^l,zeUs\^. (6.4) 



,1 k-1 

For a = 0, P{z) has the following behavior as z —>■ 1: 

^(^)=o(!:^|::l| log|^:l|)' -^-M^c/AE. (6.5) 

For q; > 0, P{z) has the following behavior as z — > 1: 

O ^j, as2;— >1 outside the lens, 
O I ! A a ^ j , as z ^ 1 inside the lens. 



To find P, we concentrate first on parts (a), (b), and (d). Later we will also consider 
the matching condition (6.3). 

Recall that h is analytic with positive real part on U. Hence, the function 

W{z) = {{z - 1)"(^ + lfh{z)) (6.7) 
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is defined and analytic for 2; G f/\ (—00, 1]. Tlie brancli of tlie square root is cliosen wliich 
is positive for z > 1. Note that 

{e°''^^w{z), forlm^>0, 
(6.8) 
e-^^iw{z), forlm^<0. 

Thus in particular, W^{x) = e^°'^'^w{x) for x e (— 1, 1), and so 

W+{x)W_{x) =w{x), for x G (-1,1). (6.9) 

Thus W satisfies the same jump relation as the Szego function D[z). Note however that 
W does not extend to an analytic function on C \ [— 1, 1]. 
We seek P in the form 

P{z) = En{z)P'-^\z)W{z)-'''ip{z)-'^''', (6.10) 

where the invertible matrix valued function En is analytic in a neighborhood of Us- By 
(6.1) and (6.10) the matrix function P'^^^ should satisfy on Usr\{T.1[J'El) the jump relation 



P^'\z) 



w{z)-W\z) 1 



= P'_>(z) 



where in the last step we used (6.8). In e^"'^'^ the + holds on H" and the — on S3. For 
x G (1 — 5, 1), we find using (6.2), (6.9), (6.10), and the fact that i^^{x)ip-{x) = 1, 

pi'\x) = pi'^(x)(^-(x)-"-iy_(x)-- "^^^^^w^xr^^ixr^ 

(1). / w{x)W_{x)-^W+{x)-^ 



= P^_>{x) 



-w{x)-^W_{x)W+{x) 



= ^^"(^) (-1 I, 

We therefore see that we must look for a matrix valued function P^^^ that satisfies the 
following RHP. 

RHP for P(i): 

(a) P^-^\z) is defined and analytic for z G Us^^ \ S for some 60 > 6. 



(b) P^^^ (z) satisfies the following jump relations: 

(^gL °) , ior ze Us n S?, 



Pl'\z) 


= P^\z) 


< Pi'\x) 


^P^}\x) 




^P^\z) 



J), forxGC/^nE", (6.11) 
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(c) For a < 0, the matrix function P^^\z) has the following behavior as z — > 1: 

/U _ 1 |"/2 U _ 1 |«/2\ 



For q; = 0, P^^\z) has the following behavior as 2; — > 1: 



For a > 0, P'^^^z) has the following behavior as 2; — > 1: 

p^^Hz) = 



,Z - \z - l|-"/2\ 

C 1 1^ _ L _ 1 > as 2; ^ 1 outside the lens, 



I |;2-l|-"/2)' as ^ ^ 1 inside the lens. 

Remark 6.1 Note that the condition (c) in the RHP for P^^^ follows from the condition 
(d) in the RHP for P, since 

P^^\z) = E,,{z)-^P{z)W{zf^ip{z)-'"'^ 

where 99 is bounded and bounded away from near z — 1, and W{z) behaves like c{z — 
1)"/^, with a positive constant c, as 2; ^ 1. 

Remark 6.2 Observe that the jump matrices for P^^^ in (6.11) are constants on the three 
parts E^, and E^. 

To construct P^^^ we first define a function g as 

g{z) = log^(^) = log {z + {z^ - 1)^/2) , for ^ e C \ (-00, 1]. (6.12) 

This function is analytic in C\(— cxo, 1] and for x e (—1, 1) we have, since ip^{x)ip_{x) = 1, 
that g+{x) = —g^{x). Therefore g"^ is also analytic across (—1,1). Since g'^(l) = 0, the 
singularity at 1 is removable, so that the function / defined by 

f{z) = g\z)/4, for z e C \ (-00, -1], (6.13) 

is analytic in C \ (— 00, — 1]. The behavior of g{z) near 1 is 

g{z) = V2{z-iy/^-—{z-lf^ + 0{{z-lf^), asz^l, (6.14) 

and the behavior of f{z) near 1 is 

m = liz - 1) - ^{z -ir + O {{z - If) , asz^l. (6.15) 

Because / is analytic on C\(— 00, —1] and /'(I) = 1/2 7^ 0, / is a one-to-one conformal 
mapping on a neighborhood of 1. A more detailed analysis of / reveals that / is one-to- 
one on the disk I2; — 1| < 2 around 1. In any case, since 0<5<r<l, we have that 
/ is one-to-one on Us- Since f{z) is real for z real, we then have that / maps Us n C± 
one-to-one onto f{Us) n C±, where C+ = {2; | Im 2; > 0} and C_ = {2; | Im2; < 0}. 

We will use the mapping ( = n? f{z) to transfer the RHP in the 2;-plane to a RHP 
in the ^-plane. Consider the following RHP in the complex ^-plane, on a contour E^ 
consisting of three infinite rays, 

7i : argC = 27r/3, 72:argC = 7r, 73 : argC = -27r/3, 

oriented as in Figure 3. 
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Figure 3: The contour S$ 



RHP for 

(a) \I/(C) is analytic for C G C \ S^, 

(b) \I/(C) satisfies the following jump relations: 



^+(0 = ^-(c)(^L J 



for C E 7i, 



vi/+(C) = ^-(C)(_°i J), force 72, 



VI/4C) = ^-(C)(g-L J), force 73. 



(c) For a < 0, the matrix function ^^(C) has the following behavior as C 



^(C) = o 



|^|a/2 |^|a/2 



as C — ^ 0. 



J^|a/2 |^|a/2y, 

For a = 0, ^E'(C) has the following behavior as C 0: 

'fogici fogicr 



^(0 = o 



as C ^ 0. 



Jog Id log Id 
For a > 0, ^^(C) has the following behavior as C — 0: 

Oyll\a/2 |^|-"/2j' asC^0in|argC|<27r/3, 

Or^l\-a/2 1^1-/2)' asC^0in27r/3< largCI < 
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Figure 4: The definition of Ei, S2 and S3 in Us 



In tlie above, we liave cliosen tlie angle 27r/3 rattier arbitrarily. It could be replaced 
by any angle in (0,7r). However, having chosen the angle, we define the contours Si fl 
and S3 n Us in Us as the preimages of the rays 71 and 73 under the mapping ^ = n^f{z), 
see Figure 4. Note that this definition does not depend on n. Then it follows immediately 
that if ^E'(C) is a solution of the RHP for \I' we can define 



(6.22) 



and P(^) will satisfy the RHP for P^^\ 

We construct a solution \I' out of (modified) Bessel functions of order a. We refer to 
Chapter 9 of [2] for definitions and properties of the modified Bessel functions and Ka, 
and the Hankel functions Ha^ and Ha \ as well as their relation with the ordinary Bessel 
functions Ja and Y^. We define \I/(C) for | a.Tg(\ < 2n/3 as 



27r<V2j/ (2^1/2) 



^^a(2C^/^) 



(6.23) 



For 27r/3 < arg(^ < vr we define it as 
iiJ«(2(-C)V2) 



iijF(2(-C)i/2) 



*(C) 



TTC 



C^/^f/jW) (2(-C)'/') (2(-C) 



r(2) 



(6.24) 



And finally for — tt < a.Tg( < — 27r/3 it is defined as 



|^^^(2(-C) 



*(C) 



l/2\ 



-nC^' [hP)' (2(-C)1/^) vrC^/^ (ifP)' (2(-C)^/^) 



(6.25) 



Theorem 6.3 The matrix valued function ^(C)? defined in (6.23)-(6.25), is a solution of 
the RHP for ^. 

Proof, (a) Since the functions la, K^, Ha \ and Ha^ (as well as their derivatives) are 
defined and analytic in the complex plane with a branch cut along the negative real axis, 
the matrix valued function \l/ defined by (6.23)-(6.25) is analytic in the respective regions. 

(b) We need to verify that the jump conditions (6.16)-(6.18) are satisfied. We first 
verify the jump condition (6.17) on the negative real axis. The Hankel functions evaluated 
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at 2(— (^)^/^ that appear in (6.24) and (6.25) are analytic across the negative (^-axis. Only 
the factor C,^/"^ appearing in the second rows of (6.24) and (6.25) has different boundary 
values on the + and — sides of the negative real axis. Taking this into account, we easily 
see that the jump condition (6.17) is satisfied. 

Next, we check (6.16). For C, G 71 wc use (6.23) to evaluate ^1/+(C) and (6.24) to 
evaluate ^_(C). From (6.24) the (1, 1) entry on the right of (6.16) is 

^i...l(^Ci)(2(_C)V2) + ^^)(2(_C)V2)). 

Using (-C)^^^ = C^''^exp(-|7ri) for C ^ 7i and formulas 9.1.3 and 9.1.4 of [2], we find 
that this is equal to 

e^°^V„(2C^/^e-^'^^). (6.26) 
If we now use the fact that, cf. formulas 9.1.35 and 9.6.3 of [2], 

Ia{z) = e^"'''Ja{ze-^'''), for - tt < arg^ < |, 

we see that (6.26) is equal to /a(2C^/^), and this is the (1, 1) entry on the left of (6.16). 
So the (1,1) entries of both sides of (6.16) agree. 

The (1,2) entry on the left of (6.16) is {i/Tr)Kai2C^^). Using formula 9.6.4 of [2], we 
rewrite Ka in terms of the Hankel function Ha as follows 

^X.(2C^/^) = ie-i'^-i/^)(2(-C)^/^). 

TT Z 

This is exactly the (1,2) entry on the right of (6.17). So the (1,2) entries of both sides 
(6.16) also agree. 

Now we know that the first rows of the left-hand side and the right-hand side of (6.16) 
are the same. To see that the second rows agree as well, we observe that in both (6.23) 
and (6.24) (and also (6.25)) the second row is equal to the derivative with respect to of 
the first row multiplied by the factor 2ni(. Since the jump matrix does not depend on (, 
the equality of the first rows yields the equality of the second rows as well. Thus (6.16) 
holds. Similar considerations show that (6.18) holds as well. 

(c) The behavior near stated in (6.19)-(6.21) follows easily from formulas 9.6.4, 9.6.6, 
9.6.7, 9.6.8, 9.6.9, 9.6.26 and 9.6.27 of [2]. □ 

Having \E' we define P^^\z) = 'if{n^f{z)) as in (6.22) and P^^) is a solution of the RHP 
for P^^\ Then we define as in (6.10) 

P(z) = (6.27) 

with En an analytic matrix valued function in Us^. We do not have En yet, but no matter 
how it is chosen, P satisfies parts (a), (b) and (d) of the RHP for P. The multiphcation on 
the left by the analytic factor En has no influence on parts (a), (b), and (d). Observe that 
for part (d) it is important that the matrices in the 0-terms in (6.4)-(6.6) are constant 
along the columns. 

We use the freedom we have in choosing En to ensure that (6.27) satisfies the matching 
condition (6.3) of the RHP for P. To that end we need the asymptotic behavior of the 
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Bessel functions at infinity. Inserting the asymptotics of the modified Bessel functions 
la and Ka as given in formulas 9.7.1-9.7.4 of [2] into (6.23), we find that uniformly as 
C — > oo in the sector | arg^l < 27r/3, 



'_^^-i/VC^/^ (1 + 0{C"')) ^C-^/^e-^C^^^ (1 + 0(C-V2))^ 
^K'/'e^^"' (1 + O(C-V^)) ^C"e-'^''' (1 + OiC'/')) 

Inserting the asymptotics of the Hankel functions given in formulas 9.2.3, 9.2.4, 9.2.13, 
and 9.2.14 of [2] into (6.24) and (6.25), we find that the asymptotics of ^'(C) in the other 
two sectors are given by exactly the same formula (6.28). Thus (6.28) holds uniformly as 
C — > oo, C e C \ E^. Taking C, — n^f{z), we then have 

*(rzV(^)) = {27:nr^" fizy^^l'^ A + ^(ti) ' + ^2nf^/^^.)a, (g_29) 

V2 l^z + 0(i) 1 + 0(1)7 ^ ' 

as n ^ oo, uniformly for z e dUs\T,. Since 2/^/^ — 9 — log0 we have exp(2n/^/^(2;)(T3) = 
ip{zY''^. Combining (6.27) and (6.29) we then get 

/9 /. 1 /l + C'(i) i + Oi^)\ 

Piz)^EJz){2T:n)-''''^ fizy^'^^l f 'l ]w(z)-^\ (6.30) 

^ ^ ^ ^ ^ ' V2 l^^ + 0(i) l + 0(^); ^ ^ ' ^ ^ 

as n — > oo, uniformly for z e ^C/^ \ E. Since we want P to match with N on ^C/^, we now 
define En as 

Eniz) = N{z)W{zr^ 7) f{zr/' {27:nr/' . (6.31) 

Then En is defined and analytic in U \ (— oo, 1]. With this En, we see from (6.30) that 

P{z)N-\z) = N{z)W{zP (^I + O (^^^^ Wizy^-^z), (6.32) 

as n — > oo, uniformly for z e ^C/^ \ E. Since W as well as all entries of N are uniformly 
bounded and uniformly bounded away from on dUs, we get from (6.32) that the matching 
condition (6.3) is satisfied. It is also clear that En is invertible, since by (6.31) En is a 
product of five matrices, all with determinant one, so that 

detEn{z)^l. (6.33) 



Now everything is fine, except for the fact that En from (6.31) is defined and analytic 
on U \ (— oo, 1], but we need it to be analytic in a full neighborhood of 1. We check 
that En has no jumps across the interval (—1, 1), and in addition that it has a removable 
singularity at 1. To that end, we first determine the behavior of W/D near the point 
z = 1. The lemma states more than what is needed for the proof that En is analytic. The 
extra information will be useful later on. 
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Lemma 6.4 For z e Us, we have 

TO^/i coefficients Cn given 6y (1.11). /n particular, we have that 

^ = 1 + 0(|.-1|V2), a..^l. 
Proof. We note that by (5.4) 

.2/ N {Z-IY{Z+If ( {Z^ - lyl"^ \ogh{x) dx 



DHz) = ^ — 7 ; ^ exp ^ / 7 . (6.34) 

Let 7 be a closed contour in U going around the interval [—1,1] in the positive direction 
and also encircling the point z e Us- After an easy residue calculation, we then find 

1 \ogh{x) dx _ loghjz) 1 r loghjC) dC 

n y_i vr^^ ~ (z^ - 1) V2 2^X(e-i)V2c-^- 



Inserting this in (6.34), we get 

DHz) - ^^^^n^M.) exp (-(z^ - 1)V^^ / J^gMCL 



Since ^^^(z) = (-2 - l)''(-z + l)^/i(-z), it then follows that 



.(.)-exp((.^-l)V^^^J-^ 



Since for z E Us, 



J_ / log MO dC _^ 



with coefficients given by formula (1.11), the lemma is proved. □ 
Proposition 6.5 The matrix valued function En{z) defined in (6.31) is analytic in U \ 

(-00,-1]. 

Proof. We first show that En has no jump across the interval (—1, 1). Let x G (—1, 1). 
Since N has the jump (5.1) and W+{x)W-{x) = w{x), we have 

AT ( MI/ ( \T3 AT ( \ I W{x)\ ( W[x) ^ 



-w[x)-^ ) \W_[x)^ 
- ^^-^""'[-W.ix)-' J 

= N4x)W^^{x) (^^^ J^. (6.35) 
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Furthermore, we have that f{x) is negative, and taking fourth roots, we get — 

if-{xy^'^, so that 



i 
-i 



^ f-{xY''\ (6.36) 



By (6.31), (6.35), and (6.36) we then have for x e (—1, 1), 

{E^)^{x) = N^{x)W^{xr^(^\ U{xr'\2nnr'' 

= {En)-ix). 
Therefore En is analytic across (—1, 1). 

Next, we have to show that En has a removable singularity at 1. Using Lemma 6.4 
and (5.5) we see that 



N{z)W{zy' = D 



oo 



' a(z)+a{z) 1 a(z)-a(z) l \ , \ \ crs 

^ '^W{z)^ 



2 2i 

a(z)-a(z)-i a(^)+a(z)-i / \ D( z) 
-2i 2 



= 0(1^ !^ ll-I/IVfi ?V ^s^^i. (6.37) 



11-1/4; -\^o 1 

Since / has a simple zero at 1 we also have 



iz-ii-i/V' 



as ^ ^ 1. (6.38) 



Thus by (6.31), (6.37), and (6.38) we have that all entries of En{z) have at most a square- 
root singularity at z = 1. Since we already know that En is analytic in a punctured 
neighborhood of 1, the isolated singularity at 1 is removable. Thus, En is analytic in 
C/\ (-00,-1]. □ 

The construction of the parametrix P near z = 1 is now completed. A similar construc- 
tion based on Bessel functions of order /3 yields a parametrix P defined in a neighborhood 
Us — {z \ \z + 1\ < 5} oi z — —1 that satisfies the following RHP. 

RHP for P: 

(a) P{z) is defined and analytic for z e Usq \ S for some So > 5. 

(b) P{z) satisfies the following jump relations on Us fl E°: 

P^z) ^ P^z) (^^^^^J^^^^_,n fj, for^eC/5n(E?UE^), (6.39) 
P+i^) - P-i^) (J^^yi "^0^ forxeC/,nE°. (6.40) 
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(c) On dUs we have, as n — > oo 

P{z)N-\z) = I + O , uniformly for ;2 e \ S. (6.41) 

(d) For /3 < 0, the matrix function P{z) has the following behavior as z — > —1: 

For /3 = 0, -P(-2) has the following behavior as z ^ —1: 

P(2;)=C> , , ' , ,1 > as z ^ -1,2; e C/^ \ E. 
^ ^ yiog |2; + 1| log 1^; + l|y ^ 

For /3 > 0, has the following behavior as z — > —1: 



O ^1 , as z — > — 1 outside the lens, 

z + lt'^ 1 
z 



O [ \^ _^ ^ ) , as 2; ^ — 1 inside the lens. 



Details of the construction of P are as follows. Instead of working with (p it is more 
convenient to introduce 

^{z)^ip{-z)^-ip{z), (6.42) 

which is defined and analytic in C \ [—1, 1]. Replacing (/? by <^ does not change the jump 
matrix in (6.39). Since ip+{x)(p^{x) = 1 for x e (—1, 1) it is obvious that we also have 

ip+{x)ip-{x) = 1, for X e (-1, 1). (6.43) 

The function 

W{z) = {{1- z)°{-l- zfh{z)y^\ ior z eU\[-l,oo), (6.44) 

is analytic in U \ [—1, 00), and satisfies 

f e-^'^'wiz), for Im^ > 0, 
I e' w[z), forlm^<0. 

Thus in particular, W^{x) — e^^^'^w{x) for x e (— 1, 1), and so 

W+{x)W-{x) = w{x), for X e (-1, 1). (6.46) 
Writing 

P{z) = Er,{z)P'^'^\z)W{z)-'''ip{z)-'"'', (6.47) 

with an invertible analytic factor En, we see, using (6.39), (6.40), (6.43), (6.45) and (6.46), 
that P^^^ should satisfy the following RHP. 
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RHP for P(i): 

(a) P^^\z) is defined and analytic for z E Uso\T, for some So > S. 

(b) P'^^^z) satisfies the following jump relations: 

1 0' 





^P^\z) 




^P^^\x) 




^P^\z) 



e 



, ioT zeUsD El, 



(^_°^ J^, forxeC/^nE^, (6.48) 



1 



(c) For /3 < 0, the matrix function P^^\z) has the following behavior as 2; — > —1: 
P^^^(z) =0 r ^ + 11/3/2), asz^-l,zeUs\E. 

For P — 0, P^^\z) has the following behavior as 2; — > —1: 



^ ' \\og\z + i\ iog|z + i|y' 



as z — > —1, z & Us\Ti. 



For /? > 0, P^^\z) has the following behavior as ^ ^ —1: 

*^V|^+l|^/2 |^ + i|-/3/2), as ^ ^ -1 outside the lens, 



O J , as 2; — > — 1 inside the lens. 

\\z + i\-i^'^ |^ + i|-^/2y' 



The jumps (6.48) of P^^^ are analogous to the jumps (6.11) of P*^^\ To construct P'^^^ we 
first define a function g as 

^(z) =log(^(z), for 2 e C\ [-1,00). (6.49) 

and a function / by 

~f{z) = ~g\z)/A, for ^ e C/A [-1, 00). (6.50) 

Then C, — f[z) is a one-to-one conformal mapping of Ug onto a neighborhood of C = 0. 
However the ordering is reversed. Thus if 2; < —1, then C > 0, while if 2; > —1, then 
C < 0. The contours Ei and E3 are chosen such that /(Ei fl Us) is part of the ray 73, and 
/(E3 n Us) is part of the ray 71 in the C-plane. 

Then we need a fixed matrix function ^! defined in the complex ^-plane that satisfies 
the following RHP on the contour E^ shown in Figure 5. 
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Figure 5: The contour 



RHP for 

(a) ^(C) is analytic for ^ G C \ S^. 

(b) \E'(C) satisfies the following jump relations: 



1 



1 



oy ' 



for C G 7i, 



for C G 72, 



for C e 73- 



(c) ^(C) has the same behavior as '^{C) near C = 0, if we replace a by (3. 

When comparing this RHP with (6.16)-(6.18) we have to realize that the directions on 
the contours are reversed. The reversal of the directions is taken care of by conjugating 
with the Pauli matrix a^. Thus we take 



^(0 =cr3^(C;«^/3)cr3, 



(6.51) 



where \E'(-; a — /5) is the matrix valued function defined by (6.23)-(6.25), but with Bessel 
functions of order (3 instead of order a. Then P^^\z) = '^(n'^f{z)) is a solution of the 
RHP for P(i) and thus 



P{z) = Er,{z)^in^fiz))^iz)-^'^'-'Wizy 



(6.52) 
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with En an analytic matrix valued function, satisfies the parts (a), (b), and (d) of the 
RHP for P. As before we construct the analytic factor E^ so that the matching condition 
(c) is satisfied as well. We define En by 

^„(z)=7V(^)#(zr-^Q jj/>r/^(27rnr/% iorzeUs. (6.53) 

Similar reasoning as before shows that En is analytic in U\[l, oo), and that the matching 
condition is satisfied. 

This completes the construction of the local parametrix P near —1. Note the slight 
difference (apart from the tildes) between En defined in (6.53) and En defined in (6.31). 

For later use, we remark that the following analogue of Lemma 6.4 holds. Its proof is 
similar. 

Lemma 6.6 For z & IJs, we have 
'W{z 



D(z) 



- (^(^)"+^ exp (z^ - 1) V2 J2dn{z + 1)" , 



n=0 

with coefficients dn given by (1.12). 

7 Third transformation S R 

In the final transformation of our original problem we define the matrix valued function 
Rhy 

R{z) = S{z)N-^{z), ioT zeC\(Usl^Usl^T.), (7.1) 
R{z) = S{z)p-\z), forzGt/5\S, (7.2) 
R{z) = S{z)p-\z), iorzeUs\^- (7.3) 

Remark 7.1 The inverses of the matrices N, P, and P used in (7.1)-(7.3) exist, since 
the determinants of these matrices are equal to 1. For A^, this is easy to see from its 
definition (5.5). To sec this for P, we note that by (6.27) and (6.33) it is sufficient to 
show that det ^(C) = 1 for C G C \ It is clear that det ^(C) is analytic in C \ H^,. 
Since the jump matrices in (6.17)-(6.18) have determinant one, it follows that det \1'(C) 
has an analytic continuation to C \ {0}. By (6.19)-(6.21), the behavior of det^(C) as 
C^Ois 

o(icr), if«<o, 

0((log|C|)2), ifa = 0, 

0(1), if a > and I arg(C)| < 27r/3, 
0(|C|-") if a > and 27r/3 < | arg(C)| < tt. 



det^(C) 



Then the singularity at is removable, which is obvious in case a < 0, and which follows 
as in the proof of Lemma 4.1 in the paragraph after formula (4.13) in case a > 0. 
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Figure 6: System of contours for R 



Thus det^'(C) is entire. From the asymptotics at infinity, given in (6.28), we see that 
det ^'(C) — > 1 as C — ^ oo. By Liouville's theorem we then find that det ^'(C) = 1 for every 
(. As aheady noted this imphes that detP(z) = 1. In the same way, it follows that 
detP(;z) = 1. 

Now R given by (7.1)-(7.3) is defined and analytic in C \ (S U dUs U dUs). Then R 
has jumps on a system of contours shown in Figure 6. However, since by construction 
the jumps of S and on the interval (—1 + 6,1 — 6) agree, it follows from (7.1) that 
R has the identity as jump matrix on this interval, and so has an analytic continuation 
to (—1 + 6,1 — 6). Also, since the jumps of S and P on Us H S° are exactly the same, 
it follows from (7.2) that R has an analytic continuation across these contours as well. 
Similarly, it follows from (7.3) that R is analytic across Us fl S°. What remains are the 
jumps on a reduced system of contours, which we will call S/j, see Figure 7, and possible 
isolated singularities at ±1. We prove that these singularities are removable. 

If a < 0, we use the fact that detP(z) = 1, and (6.4) to conclude that 

p-^(z) =0 as^^l. 
Then from (4.6) and (7.2) it follows that 



Since R is analytic in f/^yjl} and a > —1, it follows that the singularity at 1 is removable. 
If a = 0, it follows in a similar way from (4.7), (6.5), and (7.2) that 

\{\og\z - 1\) {\og\z - 1\) J 

from which it also follows that the singularity at 1 is removable. 
And finally, if a > 0, it follows from (4.8), (6.6), and (7.2) that 



R{z) 
























z-l\ 


—a 


\z — 


1 


z-l\ 


—a 


\z — 


1 



as 2; 1 outside the lens, 
O ( r V, ^ r ^! „, 1 , as z — s> 1 inside the lens. 
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Figure 7: Reduced system of contours S^^ 

From this it also follows that the singularity at 1 is removable, cf. the argument given 
after formula (4.13) in Lemma 4.1. 

So in all cases we proved that the singularity at 1 is removable. Similarly it follows 
that the singularity at —1 is removable. Hence R is analytic in C \ S^j, where 

S« = dUs U dUs U ((Si U S3) \ {Us U Us)) . 

The contour S^ is oriented as indicated in Figure 7. 

We next calculate the jumps of R. For z G S^ \ {dUs U dUs), we use (4.3) and (7.1) 
to obtain 



z 



On the boundary of Us we have R+{z) = S{z)N^^{z) and R-{z) = S{z)P~^{z) so that 

R+{z) = R_{z)P{z)N-\z), for z e dUs- 
Similarly, on the boundary of Us we have 

R+{z) = R^{z)P{z)N-\z), for z G dUs- 
So we see that i? is a solution of the following RHP: 

RHP for R: 

(a) R{z) is analytic for z G C \ S/j. 

(b) R{z) satisfies the following jump relations on S;^: 

R+{z) = R4z)P{z)N~\z), hizedUs, (7.4) 

R+{z) = R^{z)P{z)N~\z), for zedUs, (7.5) 



R^iz) = R.iz)Niz) ljN~\z), (7.6) 

for z G J:R\{dUsUdUs). 
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(c) R{z) ^ I + [- ) , asz^oo. 



Remcirk 7.2 Since log > c > on the upper and lower lips of and N(z) does 

not depend on n, the jump matrix in (7.6) tends to the identity matrix at an exponential 

rate 

R+(z) = R-(z) {I + O {e-^'^)) , as n ^ oo (7.7) 

uniformly on Ei? \ {dUs U dUs). By (6.3) and (7.4) we have on the boundary of Us as 
n ^ oo that 

R+{z) ^ R4z) (^I + O (^^^^ , as n ^ oo, (7.8) 

uniformly on dUs- At the two points where dUg intersects with the arcs Ei and E3, the 
boundary values are taken within a single connected component of C \ E/j. Since R is 
continuous up to the boundary from within each component, (7.8) holds everywhere on 
dUs. Similarly by (6.41) and (7.5) 

R+{z) ^ R4z) (^I + O (^^^^ , as n ^ 00, (7.9) 

uniformly on dUs. Prom (7.7)-(7.9) it follows as in [11] that 

R{z)^I + o(^^, as n ^ 00 uniformly on C \ Eij, (7.10) 

The relation (7.10) is the main term in the asymptotics for R and it is enough to give 
the leading term in the asymptotics of the orthogonal polynomials. In the next section 
we will give a complete asymptotic expansion for R. 

8 Asymptotics for R{z) 
8.1 Asymptotic expansion for A 

We write the jump matrix for R in the form / + A(s), so that by (7.4)-(7.6) 

P(s)N-\s), for s e dUs, 

Pis)N-\s), ioTsedUs, ^g_^^ 

{w-\syp{s)-'- 1) ' ^ ^« \ ^ 

We show that A has an asymptotic expansion in powers of 1/n of the form 

A(.) ~ £ ^ as n ^ 00, (8.2) 



I + A{s) = < 



k=l 



37 



uniformly for s G Sr. On S/j \ {dUs U dUs), the jump matrix is the identity matrix plus 
an exponentially small term, so that for every A; e N, 



Afc(s) = 0, for s e \ (dUs U dUg) . 



(8.3) 



We now determine Ak{s) for s G dUs- The modified Bessel functions Ia{C), Ka{C), 
and their derivatives, have an asymptotic expansion as ^ — >■ oo, see formulas 9.7.1-9.7.4 
of [2]. Using this expansion in (6.23), wc find that 



*(C) = (2^)~^^^'C-'^^/';^ 1 



1 i 



J y> (a,A:- 1) 

^ 4^=^/^ ^ (-l)n^-|)i i(a^ + i^-|), 

as C — oo, uniformly in the sector | argC| < 27r/3. Here (a, 0) = 1 and 

(4q;2 _ i)(4q;2 _ 9) . . . (4^2 _ - 1)2) 



^2^1/2 



0-3 



U) 



(a, k) 



see also e.g. [3, §4. 9, §4. 12]. Inserting the asymptotic expansions of the Hankel functions 
into (6.24) and (6.25), we obtain the same formula (8.4) as C — ^ cxd in the other sectors. 
Thus (8.4) is valid uniformly as ^ — > oo. Then by a calculation using (6.12), (6.13), (6.27), 
(6.31), and (8.1), we find that uniformly for s G dUs, 



A{s) = P{s)N-\s)-I 



N{s)W{s) 



0-3 



(a, k — V 



Ik _ 1^ 



j^^n^2'^g{s)'^ y (-l)^(A;-i)^ 



-{k-\)i 



Thus, 



A,(s) 



(a, k — 1) 



N{s)W{s) 



<^3 



Wisy^N-^is), for s G dUs, 
Analogously, using (6.49)-(6.53), and (8.1), we find 



-{k-\)i 



(8.5) 



Afc(s) 



(/?,A:^1) 



N{s)W{sy 



Wisy^N-^is), for s G dUs. 
To summarize, we have shown 



(8.6) 



Lemma 8.1 A(s) has an asymptotic expansion in powers ofl/n of the form (8.2), which 
is valid uniformly for s G Sr. For s G Sr \ {dUs U Sf/^), we have Afc(s) = for every k. 
For s G dUs and s G dUs, we have that Ak is given by (8.5)-(8.6), respectively. 
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For later use, we need 

Lemma 8.2 (a) For every A; e N, the restriction of to dUs has a meromorphic 
continuation to a neighborhood Usq of 1 for some Sq > S. The continuation is 
analytic, except at 1, where it has a pole of order at most [{k + l)/2]. 

(b) For every A; e N, the restriction of to dUs has a meromorphic continuation to 

a neighborhood tls^ of —1 for some Sq > 5. The continuation is analytic, except at 
—1, where it has a pole of order at most [{k + l)/2]. 

Proof, (a) The formula (8.5) shows that Aj. is analytic on U \ (— oo, 1]. We recall, see 
(6.35), that for x e (-1,1), 

N+{x)W+{xy' = N^{x)W-{xy' (^^^ . 

Since g+{x) — —g-{x), we then have for x e (— 1, 1), 



-1 0/ \ + VI 



W-ix^^NZ^x) 



1 



I iiF -T-— n — I in — p7 

k I k ^ 2 a! 2) 



W^ixy^NZ^x) 

= Afc,_(x). 

Thus Afc is analytic across (—1, 1). 

To determine the behavior near 1, we note that g{z) = 0{\z — as ^ — >• 1, and 

as z — > 1, see (6.37). Thus we have by formula (8.5), 

A/\ \z — 1\ 2 2 \ 

Ak{z) = 0\\ as 2; ^ 1. 



Since z = 1 is an isolated singularity of Ak{z), it then follows that z = 1 is a pole of order 
at most [{k + l)/2]. This proves part (a). 

Part (b) is proved similarly. □ 
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8.2 Asymptotic expansion for R 

As in Theorem 7.10 of [11] we obtain from (8.2) that R has an asymptotic expansion of 
the form 

R{z)r^I + f2^^ asn^oo, (8.7) 

k=l ^ 

which is vahd uniformly for 2; G C \ {dUs U dUs)- We have that 

Rk{z) is analytic for z e C \ {dUs U dUs), (8.8) 

and 

Rk{z) = 0[-], as^^oo. (8.9) 



The expansion (8.7) is valid uniformly near infinity in the following sense. 

Lemma 8.3 For every I there exists a constant C > 0, such that for every z with \z\ > 2, 



fc=i 



Here \\ ■ \\ denotes any matrix norm. 

Proof. Prom i?+ = i?_ + i?_A we obtain the following integral representation for R (cf. 
the discussion in the proof of Theorem 7.8 in [11]), 

R(z)^I + ^ I ^-(^)^(^) ds, forzeC\E^. (8.11) 



27ri ./sfl s- z 



Expanding the jump relation = i?_(/+ A) using (8.2) and (8.7), and collecting powers 
of 1/n, we find for k = 1, 2, . . ., 

k 

Rk,+ {s) - Rk,-{^) = Yl ^k-j,-{^)Aj{s), for sedUsU dUs, (8.12) 

where i?o,- — I- Prom the Sokhotskii-Plemelj formula, (8.3), and (8.9) it then follows 
that 

h 

Mz) = ^ / 5^^fc-.,-(s)A,(5)^, for ;2 e C \ (8.13) 

JT.R * ^ 

Combining (8.11) and (8.13) we then have for 2; e C \ E^^, 



k=l '^^R \ k=l j=l 



Rk-j,-{s)Aj{s)\ ds 



I s — z 
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One easily verifies that we can rewrite the term in the parenthesis as follows 



k=l 3=1 \ k=l 

+ (ms) - 1 ^) + E ^ i: . (8.14) 

V k=l / k=l \ j=l / 

From (8.2), (8.7), and the fact that each Aj is bounded on T,r, it follows that each term 
in (8.14) is O (n^*^'+^)) as n — > oo, uniformly for s E S/j. Hence the norm of the left-hand 

side of (8.14) is < — ^ for some constant c uniformly for s e Hr, and so 

77" 



c 1 f Ids I 

< 



n'+^ 2n ly„ \s — z\ 



Since the contour S/j is bounded and |s — -z] is uniformly bounded away from for s G 
and \z\ > 2, we obtain (8.10) for some constant C. This completes the proof of the lemma. 

□ 

The functions Rk can be computed from the known form (8.2), (8.5) and (8.6) of the 
jump matrix. We will do the calculations explicitly for the first two functions Ri{z) and 
R2{z). It will then be clear how to proceed for the further terms. 

8.2.1 Determination of Ri 

Expanding the jump relation R^ — R_(I + A) up to order 1/nwe find, cf. (8.12) 

Ri,+is) - Ri,-is) = Ai(s), for sedUsU dUs- (8.15) 

Then (8.15) and (8.8)-(8.9) with = 1, give an additive RHP for Ri. This RHP can 
be solved by the Sokhotskii-Plemelj formulas, but in this case wc can state the solution 
explicitly as follows. Prom Lemma 8.1 we know that Ai(^) is analytic in a neighborhood 
ol z — —1 and z —1 except for simple poles at ±1. So 

Ai(^) = r + 0(l), asz^l, Ai(^) = ^- + 0(1), as^^-1, 

Z — i. Z + 1 

for certain constant matrices A^^^ and B^^^ . By inspection we then see that 



AW 5(1) 
+ 



for ^ e C \ (Us U Us) , 



RJz) = { ^ + ' ' (8.16) 

- + ^— --Ai(^), for ^ e C/5 U 

z — \ z + 1 

solves the additive RHP (8.8), (8.9), (8.15). 

The constant matrix can be calculated from the formula (8.5) with k = 1, and 
the known behavior of g, N, and W near 1. The behavior of g is given by (6.14), and for 
N and W we have as in (6.37) 



a{z)+a{z) ^ 


a(z)—a(z) ^ 


2 


2i 




a{z)+a{z)~^ 


-2i 


2 




41 



N{z)W{zr = ( I ( ^ ) , (8.17) 



W{z) 



0-3 



where W/D has an expansion near 1 given in Lemma 6.4. The functions a{z) — (^^|ji/4 

and ip{z) = z + {z'^ — 1)^/^ that appear in (8.17) and Lemma 6.4 also have expansions 
near 1. We then obtain 

^(') = ^^D2(7 (8.18) 

Similarly, B^^^ follows from (8.6), Lemma 6.6, and the behavior of g, N, and W near 
— 1. The result is that 

This completes the determination of 

Remark 8.4 To calculate A^^^ we only need the leading terms in the expansions of g, N, 
and W. So in particular it is enough to use W{z)/ D{z) = 1 + 0{\z — 1|^^^) as z ^ 1, see 
Lemma 6.4. Similarly, for B^^^ it is enough to know W{z)/D{z) = 1 + 0{\z + Ij^^^) as 
z — > —1, see Lemma 6.6. This means that the analytic factor h in the weight only plays 
a role through the constant in the determination of Ri. 

8.2.2 Determination of R2 

Expanding the jump relation R^ — R_{I + Ar) up to order l/n^ we find, cf. (8.12), 

i?2,+(s)-i?2,-(s) = i?i,-(s)Ai(s)+A2(s), iorsedUsUdUs. (8.20) 

We view (8.20) as a jump relation for R2 on dUs U dUs with additive jump matrix 

i?i,_(s)Ai(s) + A2(s). 

Then (8.20) and (8.8)-(8.9) with k = 2, give an additive RHP for R2. 

Since i?i _(s) is the boundary value of the restriction of Ri to the disks Us and Us, 
and since Ai and A2 are analytic in a neighborhood of —1 and 1, except for simple poles 
at ±1, see Lemma 8.2, we have 

Ri{z)A,{z) + A2(^) = ^ + 0(1), as ^ ^ 1, 

z — 1 

and 

S(2) 

Ri{z)A,{z) + A2{z) = —— + 0(1), as ^ ^ -1, 

z + 1 

for certain constant matrices A^'^^ and B^'^\ As in the determination of we then find 



^(2) 5(2) 
+ 



for z e C \ (Us U Us) , 



— r + — r - Ri(z)&i{z) - A2(z), for zeUsU C/j. 

The calculation of A^"^^ is somewhat more complicated than the determination of A^^^ . 
It involves the two term expansion of the functions g, W/D, a, and ip near z = 1. In 
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particular we need two terms in the expansion in Lemma 6.4, which means that the 
constant Cq comes in. Using Maple we are able to compute A^'^\ and the result is 

(2) _ - 1 f A2{a,p,co) iB2{a,p,co)\ 

~ 256 ^zC2(a,/?,Co) ^^(a, /3, Cq) ' ^^"^^^ 

with Co given by formula (1.11) and 

A2{a,f3,co) = 8a + 8/3 + 8co- 4/5^ + 1, 
D2{a,p,Co) = -8q; - 8/3 - 8co - 4/?2 + 1, 
B2{a, /3, Co) = -8a - 8/3 - 8co + 4q;^ + 4/3^ - 10, 
C2(a,/3,co) = -8q; - 8/3 - 8co - 4q;2 - 4/3^ + 10. 
Similarly, the result for B^"^^ is 

r{2) _ 4/^^ ~ 1 n'^s /^-^2(/3, a, rfo) iS2(/3, a, do) \ r.-'^s (o 9Q^ 

~ 256 l^iC2(/3,Q;,do) -D2{P,a,do)J ' 

with do given by formula (1-12). 

Remark 8.5 To calculate R2 we need two terms in the expansion of W/D near 1 and 
W/D near —1, cf. Lemmas 6.4 and 6.6. That's why we see the influence of the analytic 
factor h in the form of the coefficients cq and do- 

8.2.3 Result for general Rk 

In a similar (although increasingly more complicated) way we can determine Rk for every 
k. The calculations become unreasonable to do by hand. Assisted by the computer 
algebra package Maple, we find for i?3, 

^(3) 5(3) (J{3) 2^(3) _^ 

^'W = I3T + lTT+(^+(ITTF- fo-ec\(c/.uc/.). (8.24) 

where ^1^^^ and B^^^ are given by 

^(3) ^ l^!^ll^<^3 ( ^3(tt, /3, Co, do) iBs{a, /3, cq, do)\ j^-^s 



8192 °° \^iC3(Q;,/3,co,(io) D3{a, (3,co,do) 



oo 



8192 °° V^<^3(/3,a,do,co) -D3(/3, do, cq)^ 
with Co and do given by the formulas (1.11) and (1.12) respectively, and 

As(a,(3,co,do) = 16(4/3' - l)(co + do + 2q; + 2/3) - 40;' (4/3' - 1) 

-128(q;' + /3' + cl) - 256co(q; + /3) - 256q;/3 + 12/3' - 16/3^ - 2, 

L'3(a,/3,co,do) = 16(4/3' - l)(co + do + 2a + 2/3) + 4ai' (4/3' - 1) 

+128(q;' + /3' + eg) + 256co(q; + /3) + 256q;/3 - 12/3' + 16/3^ + 2, 
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and 

B3{a,P,co,do) = 384(q; + /3)^ + 384co(co + 2q; + 2/3) - 388q;^- 252/3^ + 64q;^ 

+48/3^ + lUa^(3^ + 528 - 384(« + /3)(a^ + l3^) + 960(a + l3) 
-192/3^(00 + do) - 384coa^ + 912co + 48do, 

C3{a,P,Co,do) = 384(q; + /3)2 + 384co(co + 2q; + 2/3) - 388q;^- 252/3^ + 640;^ 

+48/3^ + lUa'^p^ + 528 + 384(a + /3)(a^ + (3'^) - 960(a + /3) 

+192/3^(co + do) + 384coa^ - 912co - 484- 

For C(=^) and D(-^\ we find 

(3) ^ (4a^-l)(4a^-9)(4a^-25) /-I z 
12288 V i 1 

£,(3) , 



(4/32 - l)(4/32 - 9) (4/32 - 25) f-l -i 



12288 V -i 1 

For general k, we get that it!fe(-z) in the region C\ 5 U f/^ j is a rational function with 

poles at 1 and —1. The residues in 1 and —1 are denoted by A^'^^ and B^'^^ respectively. 
Then 



z+1 \Z' 



R^(z) = - — - + + { — ) , as z^oo. (8.25) 



9 Proof of Theorems 1.4, 1.6, and 1.10 

In this section we go back to the orthogonal polynomials and prove the Theorems 1.4, 
1.6, and 1.10. By tracing back the steps 

Y ^ S ^ R 

and using the asymptotic expansion for R obtained in the previous section, we obtain an 
asymptotic expansion for Y. We will restrict ourselves to the asymptotics in the exterior 
region C\ [—1, 1]. The asymptotics on the interval (—1, 1) and near the endpoints will be 
considered in the next section. Let i^T be a closed set in C \ [—1, 1]. Then we can choose 
the lens and the regions Us, Us disjoint from K. From (3.7), (4.2), and (7.1), we then 
obtain 

Y{z) =2-'"'^R{z)N{z)ip{z)'"'^, iorzeK. (9.1) 
9.1 Asymptotics for 7r„ 

We start with the proof of Theorem 1.4 concerning the asymptotics of the monic orthog- 
onal polynomials nn{z) for z E K. 
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Proof of Theorem 1.4. Since 7r„ = Yu we obtain from (9.1) 

TTniz) = (Ru(z)Nu(z) + Rl2(z)N2l(z)) . 

If we use (5.5) and (8.7), we get uniformly for z in K, 



TT^n{z) ^ Dqo a{z) + a{z)-^ 
(^(^)" ^ D(z) 2 



1+E 



{Rk)n{z) 



+ 



k=l 

a{z) — a{z)~^ 



Since 



D^D{z) -2i 

2 V2(^2_ 1)1/4' 

we obtain Theorem 1.4 with 

i a(z) — a(z)~^ 



E 

.fe=i 



{Rk)l2{z) 



a{z) + a{z) ^ 



Uk{z) = (Rk)n(z) + 



Dla(z) 



a[z] 



-,{Rk)l2{z). 



(9.2) 



(9.3) 



Using (8.16), (8.18), and (8.19) and simpUfying the expressions, we see that Hi is given 
by (1.13). From (8.21), (8.22), and (8.23), we obtain the formula (1.14) for Ha. □ 

Note that the factor in the formula (9.3) is cancelled by the factor D"^ which is 
included in the (1,2) entry of Rk (see formulas (8.16), (8.18), (8.19) for the case k = 1), 
so that Doo does not appear in the final expression for Uk- This is due to the fact 
that we are dealing with monic polynomials in Theorem 1.4. The functions Ilk{z) are 
rational functions of (f{z) with coefficients that are polynomial expressions in a, /3, and 
the coefficients Ck and dk- 



9.2 Asymptotics for jn 

We next consider the asymptotic behavior of the leading coefficient 7„ of the orthonormal 
polynomial. 

Proof of Theorem 1.6. To emphasize the n-dependence, let us use Y^'^^ to denote the 
solution (2.5) of the RHP for Y. Then we see from (2.5) that 

yir'H^) = -27ri7X(^)- 
Since 7r„ is a monic polynomial of degree n, we therefore have that 



In 



lim z-^'Y, 



27ri 2-^-00 



21 



z . 



This expression appeared earlier in [12, Theorem 3.1] and [10, Section 3.2]. Using (9.1) 
(with n replaced by n + 1) we find that for z in the exterior region. 



7^ = 



: lim 

27ri z^oo 



■,n+l 



ip{z) 



n+1 



[zRt^''\z)N^,{z) + zR^^■^'\z)N,^{z)) 



, (9.4) 
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where denotes the R matrix corresponding to n + 1. Recall that does not depend 

on n. 

We note that 

^ = 2 + o(^^^, as^^oo, (9.5) 

and 

iVn(^) = 1 + , i?5+^^(;2) = 1 + Q , asz^oo. (9.6) 
Prom (5.5) it follows that 

and from (8.7) 

?(n+l) / iRk)2l{z) 



k=l ^ ' 



as n — > oo. 



(9.7) 



Using this and Lemma 8.3 and (8.25) we find 

lim ^i^2i(^) - V f , , asn^oo. (9.8) 
^-.oo (n + 1)*^ 

Combining (9.4)-(9.8), we obtain that 7^ has an asymptotic expansion 

Taking square roots, and rearranging terms we obtain the asymptotic expansion (1-15) 
for 7„. Using the exact expressions for A^^\ B^^\ A^'^\ and B^"^^ given in (8.18), (8.19), 
(8.22), and (8.23), we obtain the explicit expressions (1.16) and (1.17) for Fi and T2. □ 



9.3 Asymptotics for the recurrence coefficients 

Recall the three term recurrence relation 

TTn+iiz) ^{z- bn)nn{z) - alnn-i{z) (9.1) 

satisfied by the monic orthogonal polynomials. In this section we give the proof of Theo- 
rem 1.10 concerning the asymptotic behavior for the recurrence coefficients a„ and 6„. 

Proof of Theorem 1.10. First, we determine formulas for a„ and 6„ in terms of Y. We 
use an argument taken from [10, Section 3.3]. Let y^"^ and be the solutions of the 

RHP for Y corresponding to n and n + 1, respectively. Since and Y^'^^ are analytic 

in C \ [—1,1] and satisfy the same jump relation on the interval (—1, 1), it follows that 

ff{n) ._ y(n+l) ^y{n)^-l ^g ^) 
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is analytic in C\ { — 1, 1}. In the same way as was done in the proof of Lemma 2.3, we can 
prove that H^'^'> has removable singularities at the points ±1. Therefore H^'^^ is entire. 
Prom the asymptotic behavior (2.2), the definition (9.2), and the fact that detF^"^ = 1, 
it follows that 

z + 0{l) 0(1 
0(1) Oiz-') 



Using Liouville's theorem we then find that 



i\ \ , as z — > oo. 



H^^\z) ^{^ ^'''^ (9.3) 
for certain constants r^, s„, and Thus 

y(n+l)^^) ^ - Sn^ yW(^). (9.4) 

The (1, 1) entry of this matrix relation gives 

-Kn+i^z) = {z- rn)7rn{z) - 27ri7^_iS„7r„_i (2;) , 

which is the recurrence relation for the orthogonal polynomials. Comparing with (9.1), 
we see that 6„ = r„ and — 2mj^_iSn, The (2, 1) entry of (9.4) shows that —2TTij^ = tn- 
Thus 

bn^Vn, al^-Sjn-l. (9.5) 

Using (9.2), (9.3), (9.5) and the fact that dety(")(^) = 1, we find that 

6. = . - Ylr'\z)Y^^\z) + Y}r'\z)Yi^\z), 

and 

al=[Ylr'\z)YS\z)-Ylr'\z)Yl^\z)) (4"H-)4"-'^(-) - 4"^(^)4""'^(^ 

for every 2; G C \ [—1, 1]. From the asymptotic behavior (2.2) of Y we then find, upon 
letting ^ — >• 00, 

6„ = hm (z - Yl^-''\z)Yj^\z)) , (9.6) 

and 

al^]imz'Yl^\z)Yi^\z). (9.7) 



z— >oo 



We use (9.6) and (9.7) to express the recurrence coefficients in terms of R. From (9.1) 
and (9.7), we get for the recurrence coefficient 

al = Jim (^R^^{z)Nu{z) + rS{z)N,,{z)) (^Ril\z)N,,{z) + R^^j {z)N,,{z)) . 
From the behavior of i?'-"^ and N at infinity (in particular (5.5)) one then easily sees that 



Z—^00 



00 



< = hm ( zR^Az) + ^ ) ( + (^) ) • (9.8) 
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For bn we find in the same way 



= lim ( 

2— »0O \ 



z — zR 



(n+1), 
11 



{z)R^i\z)). 



(9.9) 



In the asymptotic expansion for R, we may take the hmit z — > oo, by Lemma 8.3. We 



find 



^21 + -^21 



2i 



'12 



fe=l 



(9.10) 



and 



lim 



Z — Z 



Um —2; 

2^00 



(i^fc)ii(-g) 
(n + l)'^ 

(i?.)ll(^) , ^ {Rl)22{z) 



{Rl)22{z) 



1=1 



= -E 



^ (n+1) 

ik=l ^ ' 1=1 

^11 "T -Oil . "T -D' 



jfc=i 



(n + 1)^ 



+ 



22 



(9.11) 



Expanding (9.10), we find a complete asymptotic expansion for in powers of 1/n. The 
constant term is 1/4, and the term with 1/n has coefficient 



which by the formulas (8.18) and (8.19) for A^^") and 5*^^) is equal to zero. Thus 

The full expansion of a„ follows from (9.10) and gives us the expansion (1.25) with ex- 
plicitly computable coefficients Aj^. 

For bn, we find from (9.11) that the 1/n term in the expansion has the coefficient 

which is zero by (8.18) and (8.19). Thus 



O 



n^ 



The full expansion (1.26) of bn follows by rearranging terms in (9.11). 
Remark 9.1 Note that (9.7) can also be obtained directly from (2.5), since 



□ 



y2i(^) = -27rz7^_i^"-Ml + 
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and 



Then (9.7) follows easily, because 

7n-l 



In 



The formula (9.7) is due to Deift et al. [12, Theorem 3.1], see also [10, Section 3.2]. 
However the formula (9.6) is different from the one given in these references, which is due 
to the fact that the orthogonality weight is varying with n in [12]. 



10 Proof of Theorems 1.12 and 1.13 

In this section we determine the asymptotics for the monic orthogonal polynomials 7r„ on 
the interval of orthogonality [—1, 1]. We use the behavior of the first column of Y (which 
contains the orthogonal polynomials) on [—1,1], which has already been determined, see 
[29]. In [29] this behavior has been used to derive universal behavior for eigenvalue 
correlations, associated to the modified Jacobi unitary ensemble, in the bulk and at the 
hard edge of the spectrum. 

The behavior of the first column of Y in the bulk of the interval [—1, 1] is given by 
[29, Lemma 3.1]. Namely, 



Yii{x) 



( a+{x)+a+{x) 1 a+{x)-a+{x) ^ \ . , 

2 2i \ i(n arccosi;+V'(a;))o-3 [-'-) nni9l 

a+(x)-a+{x)-^ a+{x)+a+{x)-^ I ^ ^^^"^^^ 

— 2i 2 / 

for X e (—1 + 5, 1 — 5), where ip is given by (1.34). 

Proof of Theorem 1.12. Let K ho. a. compact subset of (—1,1). We can choose the 
radius 5 of the disks Us and so small that K C (—1 + 5, 1 — 5). From (10.12) and the 
facts that 

a{z) + a{z)-' ^{zf' and »(^) " _ j ^^^^'^ no 13) 

2 V2(^2_ 1)1/4' % V2{z^-iy/'^' ^ ■ ^ 

we obtain 

Yu{x) \ _ 1 6 4 



Y2i{x)J 72(1- a;2) 1/4' 



-'''''R{x)D^ 
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for X & K. Since (p+{x) — exp(i arccosx) we then find after an easy calculation 



2-'"'''R{x)D, 



0-3 
00 



cos ((n + I) arccosx + ib(x) — 7) 
X I , \ \ 1 • (10.14) 

—i cos ((n — ^) arccos x + ip{x) — j) 

Since Yn = 7r„, and since Ru{x) — 1 + 0{l/n) and R\'2.{x) — 0{l/n) as n — > 00 with 
error terms that hold uniformly for x & K, and that have a full asymptotic expansion in 
powers of 1/n, the theorem is then proven. □ 

Next, the behavior of the first column of Y near the endpoint 1 is given by [29, Lemma 
3.3]. Namely, 



Yuix)\ 1 2nn ^-n.sM^(^) T ^(narccosa;) , 



121(2^) J \J w{x) 1 1 arccos a;J^(n arccos x) J ' 

for X e (1 — 5, 1), with M{z) given by 

M{z) = Riz)N{z)W{zr^ (^^^ fizr^\ (lO.ie) 

Proof of Theorem 1.13. Wc need to determine a convenient expression for M+(x) and 
insert this into (10.15). From (5.5), (10.13) and (10.16) we have 

^+(^) = 2(i-x2)V4^^^)^- [-ic^4x)-y' ^lixy/' 



X 



It may be shown that D+{x) — ^/w{x) exp{—i^/J{x)), cf. [29, Remark 3.2]. We then have 
by (6.8) 



W+{x) 



D^x) 

Since (p+{x) = exp(i arccos x) we then obtain from (1.36) and (10.17) 

- 2(l-.^)V4 ^(")^g (-^e<^(^) e-<^(^) ) [-^ l) ^^^^^ 
D^^^na3/^ cosCi(a;) sinCi(a;) ^ ^ ^ xa3/4 



(1-^2)1/4 °° \^-icosC2(a;) -isinC2(a:;) 

Inserting the fact that 

g7ri/4 

/+(a;)^^^ = — ^(arccos a;) "^/^ 
v2 
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into (10.18) we find 

Ja{n arccos x) 



M+{x) 



I arccos xJ'^ {n arccos x) 



^ (arccos a;)V^ ^r.-^,, f cosCi(x) sinCi(a;) Wj«(n arccos a;) \ 

^J2{1 — x'^y/'^ °° \—i cos C,2{,x) —ism(2{x)J \J'a{n arccos x) J ' ^ ' ' 

If we now insert (10.19) into (10.15), and use tlie facts that Yn = 7r„, Rii{x) = l + 0(l/n) 
and Ri2{x) = 0{l/n) as n — >• oo, we obtain (1.35). Note that the error terms hold 
uniformly for x & (1 — 5, 1) and have a full asymptotic expansion in powers of 1/n. This 
proves the theorem. □ 
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